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I. INTRODUCTION 

Tie ImStallatiOoll CGE 7a5 tuüLDIDe enjines X9 a Snif raises 
several prcblem areas in the design or the intake and 
exhaust ducting. The problems relate mainly witn the large 
volume of ccmbustion air required and the properties of the 
exhaust gases rejected to the atmospnere at іда tempera- 
tures and velocity. Вот сошрагъзоп, а borler's compustion 
air requirement is nearly stoichiometric but the gas turbine 
operates at about 400 percent of stoichiometric. Tie toil- 
er's exhaust is about 400 degrees F arter leaving the last 
rows of the economizer, but gas turbine exhaust temperatures 
are frequently as high as 950 degrees F. 

In addition to the air that passes tnrough tne gas 
turbine engine there is also a reguirenment to ventilate the 
enLgine enclosure. ÀAn adequate and uniformly distrikutedi 
cooling airflow is required around tne engine to maintain 
engine-uountel components at their proper operating tempera- 
tures ard to minimize the heat rejected to the engine room 
thereby reducing the heat exposure of perating personnel. 
Many current designs branch the engine cooling airflow off 
the main intakes and/or join neateji enclosure cooling air 
irto the engine exhaust ducting. Figure 1.1 shows a typical 
Vout of inlet and exhaust ductiag. Since the enclosure 
EsoPing airflow is on the order of 20 percent ог the 
engine's full power airflow rate, it is an important part of 
the ducting design. 

Ihe fundamental requirement of an intake design is to 
provide air to the engine compressor witi tne finimua total 
pressure loss and witn a minimum of total pressure distor- 
Eon. The loss of total pressure in the intakes leads to a 


loss of engine power and an increase in svecific fuel 
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Figure 1.1 Typical Shipboard Inlet and Exhaust Ducting. 


consumption. Schwieger reports "Typical exchanje rates are 
that a one percent loss in intake pressure is e,juivalent to 
a 2.2 percent loss in power and a 1.2 percent increase in 
specific fuel consumption" (Ref. T] Aaĉitiənally, ТӨСІ 
pressure distortion at the compressor face can lead toa 
risk cf compressor blade failure. 

Exhaust ducts must also operate ида minimum pressure 
loss. "The exchange rate is 1.1 percent loss in power and 
1.1 percent increase in specific fuel consumption for the 
one percent increase in total pressure at the power turbine 
exit" [ Ref. 1]. 

Conflicting with the desigm objective to reduce 655526 


in the ducting system are several possible requirements to 
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install components in the ducting system which contribute to 


the losses kut not directly to engine performance. Filters 
are installed to increase engine life. Silencers are 
installed to reduce noise. Macalnery arrangements dictate 


the use cf certain elbows, contractions, and transitions. 
The intared signature of the ship's exhaust plume can be 
reduced by the installation or an eductor system at the 
exhaust exit. Ihe eductor also improves the environment oc 
nast mounted equipmert and may contribute to flight safety 
when operating helicopters. Some systems use an eductor 
arrangement installed at the exhaust plane of the engine to 
pump cooling air through the engine enclosure. A waste neat 
recovery boiler may te installed in the exhaust to improve 
overall efficiency. То reduce pressure losses every attempt 
Should be made to reduce the velocit; in the duct. Lower 
melocities recuires larger ducts. Бас о: Ее” за ав 
conpromise must balance the large volume of the ship occu- 
pied with inlet and exhaust ducts and the volume for other 
uses such aS weapons and habitibiiity. In Summary there are 
many different components that can be utilized within the 
ducting system and have various effects оп the system 
performance. The effects also vary with the operating point 
ot the systen. 

It is not a straight rorwdrd problem to predict how 
components in the ducting system will perform. It is an 
Enteracting or matching type of probleu. Furthermore, it is 
a dynamic problem as parameters affecting сеггогпапсе саг 
Vary over a wide range. For example, one power setting of 
the gas turbine requires a different mass flow rate of air 
tnan  ancther. The variable mass flow rate tnrough the 
ducting system creats a variable inlet and exbaust  Guct 
pressure loss. The variation in exhaust temperature arfects 
the losses in the exhaust duct. Ustimateiy all losses 


affect the performance of the gas turbine engine. 
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One approach to the analysis of ducting system performs 
ance is to separate the problem into two areas of concern. 
The first area should deal with a one-dimensional analysis 
of the ducting system to determine now pressure losses 
affect engine performance and now the various components of 
the system contribute to the sum totai of these losses. Тһе 
second area should deal with the distortion of total pres- 
sure acrcss any section of the duct. This area becomes a 
three-dimensional  prcblem where interest is directed to 
perfcrmance not just at any section of the duct but to 
within that section to the variation of velocity across tne 
Cutting plane. Tne one-dimensional and tnree-dinmensioral 
areas of the analysis are of course related. 

The relationship between the one-dimensional and the 
three-dimensional aspect of the problem is understood and is 
dealt with in an empirical manner. Іле method is to apply a 
соггесетоп factor to the loss developed in the one- 
dimensional analysis of a particular system component, based 
on the distortion of the rlow assumed to re presented to tne 
component. ТЕ the assumptions about flow distcrtion are 
made and are accurate much valuable information results fron 
the one-dinensional analysis. 

The three-dimensional analysis of a duct system is 
possible only ‘for a very simple system and requires very 
large computer assets. It is current practice to dea xm 
thoee-dimensional analysis ort complex systems through model 
Studies. One-dimensicnal analysis on the other hand is well 
Suited for analysis cr a computer. 

It is the intent of this study to develop the nethod- 
ology for a one-dimensional analysis of a Jas turbine 
engine's inlet and exhaust ducting as mignt be instailed on 
a ship. Then to impiement the method in an interactive 
computer program which allows rapid input of the duct geom- 


веску; desired operating point and ambient conditions to 
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obtain an accurate estimate of performance. The designer 
can then decide to make changes to components to achieve 
design objectives and паке tnose chanjes to the duct yeon- 
etry tnrougn an editing routine and rerun the problem. Once 
tre designer is satisfied with the one dimersional analysis 


a firm oasis exists to provide a design for model studies. 
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II. THEORY AND ANALYSIS 


А. GENERAL 


А опе dimensional analysis of the riow in duct sections 
deml izes the Sernoulli Fquation modified to account for 
losses. The term one-dimensional 1i3 an adjective often 
applied to flow situations. The whole flow is considered to 
Ее one large streamtube with average velocity У at each 
Cross section. Thus the one dimension is the location down 
the duct. Losses refers to the pressure loss caused by 
frictional stresses in the airfiow boundary layer and by 
turbulence. A thorcugh understandiny of these terms and 
concerts is required to convey the meaning o£ the results of 


the duct system analysis. 


В. THE BERNOULLI EQUATION 


The Pernoulli Eguation is discussed in any basic text on 
fluid mecharics. It was developed to describe tne flow work 
of an ideal incompressible fluid in steady flow through a 
streantuEe. In words it states that the mechanical energy 
per unit mass along a streamline 15 conserved. Тһе 


EErproulli cyuation is: 
ue/2g. + P/p TT (G70) 2 = constant.: (egn 2:1) 


It relates velocity, pressure, and potential height. The 
constant may have a different value for each streamline, but 
for the purposes of duct flow certain simplifying assıap” 
tions are valid which паке the constant valid Гог any 
Streamline. The assumptions are that the static pressure is 


constant at any point in a cross section of tne duct. The 
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next assumption is that because the system uses gases, the 
effect of variation in potential height at a duct section is 
so small relative to the other teras that its effect is 
neglected. This assumption is extenied further to include 
the change in elevation effect at any section relative to 
any otner section. 

Alternate forms of the Bernoulli Equation are obtained 
Бу multiplying through Ey either gc/3 orp. Or interest to 
gas flow and duct design is tne form optained by multiplying 
through Ly p. Applying the above assumptions the resulting 


equation is: 
рт2/29. +p = constant (еап 227 


In this form the constant .has units of fFoot-pound гове- 
feet3 and expresses the energy per unit volume flow rate. 
It reduces to pound force/feet*® or pressure. Each term in 
the expression is given a name. The velocity term is the 
velocity presure, p is tne static pressure, and the constant 
is the total pressure. In words, the total pressure at a 
point is the sun of tne velocity pressure and the static 


pressure. 


C. MODIFIED BERNOULII EQUATION 


Altaough eguation 2.2 was derived for flow along a 
streamtute of an ideal frictionless flow it can be extended 
to analyze flow througn ducts in real systems by applying 
the First Law of Thermodynamics. A good development of the 
application of the First Law of Thermodynamics to pipe flow 
is found in {Ref. 2]. It results in tne modified Bernouili 
Equation (2.3). Equation (2.3) incorporates all thé as 
tions so far and includes the term Ap. The tlow resistan 


in a system with a real fluid between stations 1 and 2 15 


represented by the total pressure loss,4p,. 
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аа Зе Ұр. - 22229. + 0р, + ар, (eqn 2.3) 


The velocity used in the moditied Bernoulli Eguation 
will be taken as the mean velocity and then this eguation 
will be assumed valid for any streamiine in the duct. 
Analytically this is not correct because there is a varia- 
tion of velocity at a duct section iron tne walls to the 
center of the duct. Ihe error introduced by this assumption 
is offset by two circumstances. Pirot vaiti turbulent- flow 
the velocity profile is nearly uniform wnich makes the mean 
velocity a good approximation of the velocity at any point 
in the cross section. Second, experimentally determinei 
loss coefficients are utilized in computations and this 
coerficient is applied using the mean velocity. Then ir the 
velocity profile іп the system matches the profile of the 
experiment, the loss will be correctly computed using the 
mean velocity. 

The computer program uses the mean velocity and computes 
it basea on nass flow rates. The mean velocity 1S computec 
from the mass flow through a sectional area and the density 
ог the zluid at the section using equatior 2.43. Density is 
computed cy the perfect jas law equation (2.5) and isa 
function oz the absolute temperature of tne gas and the 


Static pressure of the gas. 


W 





VERS » РА (egn 2.4) 
БЕРЕКЕ (egn 273) 
where p = static pressure 

Egg secohnstant 
T = absolute temperature 
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р. PRESSURE LOSSES 


There are two types of fluid losses in the "auem 
system, frictional and dynamic losses. Frictional losses 


cccur along the walls of the entire duct length and are due 


to-r uid visc5S ty. Dynamic losses result from disturb 
tne flow such as a change of direction, contraction 9% 
expansion. 

The Darcy-Weisbach equation (2.6) calculates the fric- 


tion loss for straight duc tS: 


yè 
Darcy-Wiesbach equation Ap, = f (L/D) 5 (eqn 2.6) 
с 


where Ap, = frictions ross 
ln terms of total d pressume 
= friction Басов 
т. = аск ел түп 
D = duct diameter or 


equivalent hydraulic diaaeter 
ру” 
2.9‹ 


The friction factor, f, used in computing duct losses is 





- velocity pressure 


taken frcm a correlation Гу Swamee and Jain presented ir 
[Re£. 2]. 





геа E E 
10° 25 < 10 
0. 
7 Е е — ас 8 (еап 2.7) 
СЄ 5 + ES | 5000s Res /0 


The absolute roughness factor, e, 15 taken to be 0.00015 
feet for all air duct components. For rectangular strangi 
duct sections the equivalent hydraulic diameter, Dp, is 
calculated by equation (2.8) presentei in [Ref. 3]. 
Equations 2.6, 2.7, and 2.8 are utilized in the programm ы 
computing friction losses in the stralgnt sections of the 
duet. 
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0.625 
Do - 120 cum (ech 2-8) 


(a > 529-250 


Friction losses оСсСПІ in all 216621635 not just in 
cera ight duct. There are two tecnnisues t5 arrive at the 
friction losses in these other fittings. Тпе decision about 
which tecnnique to use depends on the wnether the fitting is 
ЕШӘШІЕС or long. шата еса ыа fri-tion 15 accounted for 
by measuring the connecting sections of straignt duct to the 
center of the fitting. No attempt is made to include fric- 
mom. jn the calculation of fluid resistance» for a short 
Meetings Elbows are skort fittings. For long fittings such 
as diffusers and contractıons, friction is inciuded in tae 
computation of the flow resistance coefricient. Therefore, 
a connecting straight duct length should be measured to tne 
center of an elbow cr to the start or eni of a diffuser or 
contraction. | 

Dynamic losses are sometimes called local ог minor 
losses. In piping systems, losses due to the local distur- 


tances of the flow are often called minor losses. In very 
iong piping systems these losses are usually insıgnificant 
шиисопрагазоп with the friction іп іле lerjth considered. 
In the duct used for a gas turbine instaliation  tnese 
so-called ninor iosses actually become major losses because 
or the snort lengths usually encountered. Experimental 
results are almost always used to account for pressure 
losses through the duct fittings. sucLosuan:ormndtriron- rris 


ually civen in the form of equation 2.9. 
Ар. - Крч2/24с (egn 2.9) 
The coerrílcient K is given for tne fitting in numerous hand- 


books. Figure 2.1 shows some typical representations of the 


ШИ ornation available. 
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K=0.35 К=0.10 


Сз 
К=0.20 | K=0.40 ` 
ne EORR were reproduced fron 
ef 
eee 


Figure 2.1 Typical K Values for Fittings. 


One of the purposes of the program is to provide K coef- 
ficients for various fittings selected to cepresent duct 
compcnents. K values can vary with tne geometry of a 
fitting: For example, a long smooth radius rectangular 
elbow has a lower K value than a snort smooth radius rectan- 
gular elkow. The prcgram taxes this into account and is the 
reascn Гог the various questions about a rittinj's geometry 
in the area of the program where the user iS inputing the 
duct systen. 

Two fittings in the prograu's menu do not require gecue 
etry inputs to obtain resistance information. The “ама 


fittings are filters and the gas turkine nodule. The reason 


for the lack of questions is that tne losses are  Eased on 
manufacturer's data. Filter aanuracturerzs provide pressure 
loss data based on face velocity and tae module is based on 
the mass flow rate of cooling air. A power curve fits the 
data and the program uses the curve to aodel pressure losses 
for these fittings. 

Table I sumarizes the fittings available from the 
program's menu. Ihe fluid resistance coefficients аге 
computed by the program upon input or the required gecuetry 
Pactors for the fitting. НЕШЕ OE the duct fittings is 
accomplished interactively. The source of the model for 
eden fitting іс noted in the projram listing іп tae title 
ШЕСЕ ог the fitting subroutine. The program subroutines 
men throwgn ГІТ29 соггеѕ рога to the fittings listed it 
Fable I . à sketch of each fittinj is provided in the 
user's manual ror the progran. тре user's manual is 


Appendix C. 


Е. GAS IURBINE/SYSTEM INTERFACE 


General Electric Company, the manufacturer of the 142500 
marine gas turbine, publishes performance data for its 
engine unđer variable operatiryg conditions. (Ref. 4]. TE 
is important to understand how the shipboard engine is oper- 
ated under variable operating conditions such as duct losses 
and ambient temperature , pressure and humiljity so that the 
proper corrections may be applied to the engine performance 


parameters for these variables. 
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09 
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11 
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14 
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17 
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22 
23 


24 


TABLE I 


Fittings Available From Program Menu 


Description 

Intake shaft, rectangular 2г055 
section, side oririces, with ог 
without louvers 


Straight duct, omg 
rectangular 


Smooth radius round elbow 


Round 90 degree sejmented elbow 
with 3,4, or 5 ээ есе» 


Mitered round elbow with or 
Without concentric vanes 


Mitered rectangular elbow 


Smooth radius rectangular 
elbow 


Smooth radius rectanjular 
elbow with splitters : 


Mitered rectangKlar elbow 
with vanes 


Rectangular elbow with 
converging OF áivor ING DE 


90 degree rectangular elbows 
in a 2-shaped corrcigurpatqon 


20 degree rectanguiar elbows 
in different planes 


Eranch section of a civerging 
чуе 


Main section of a diverging wye 


Branch section o; a convergent 
муе 


Main section of a convergent wye 
Conical round diimaiser 

Piane in-line dirfuser 
Pyraaidial in-line diffuser 
Transitional diffuser 

Round contraction 

Rectangular сопегас Шеш 

Screen obstruction in duct 
Louver entrance 


Continued next page 
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25 Filter element 


26 Multi-baifie type silencer 
217 Gas turbine modules enclosure 
28 Waste heat recovery boiler 
29 Abrupt exit 


Cad 
© 


Fittümg rot listed 


Toae Shipboard operator's point of view the engine 
Should drive the skip at the desired speed whetner it is a 
hot day cr a cold day, ог if the inlet duct losses are four 
inches of water or eight. The engine is operating а1тЕег- 
ently under such conditions to produce the same horsepower 
and speed. The proper correction  ractor set to be applied 
to the tabulated data is the set for constant speed and 
horsepower. The corrections are applied in the program with 
each iteration ‘of the duct system performance calculations 
using the current values of the inlet and exhaust duct 
losses and ambient conditions. The correctionS are very 
smail (less than twc percent) and tne convergence of the 
correct engine operating point and duct losses created by 
tne mass flow of air required at the operating point is 


guite stable. 


Е. FAN/SSYSTEM INTERFACE 


Ше орегавт по рсіпЕ ОП the ran installed іп а duct 
system is the point where the fan characteristic curve 
intersects the system characteristic curve. The fan curve 
Shows pressure rise vs. flow rate.  wWitn increasing flow the 
pressure rise across the fan is reduced. ne system curve 
is the opposite, increasing flow in tne system increases the 
resistance to  £low. Figure 2.2 represents tais situatior 


Em3phically. 
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In the iteration process the system curve js estim 
as a quadratic fitted to the origin as a  mipimua por m г 
the other point at the assumed flow and tae resulting pres- 
sure loss. Similiarily the fan curve is also represented as 
a quadratic with a maximum ас пахашип pressure attainacle 
and the corresponding flow and another point at zerom m и 
sure and maximum flow. The representation of the far 
performance for the default condition, the Spruance class 
destroyer module cooling fan, is excellent. 4ith an equa- 
tion for botna curves the point of intersection "cam 
obtained. The resulting flow is used in the next iteratioz 
until the resistance of the system and tne pressure rise 


across the fan is the same for the assumed flow. 


с. JUNCTIONS OR WYES 


An excellent discussion of tne @1x1ng of two streams 
moving at different velocities was written by Idel'chik and 
is presented here to develop the background for the 
eductor/system interface discussion. 


_ The junction of two parallel streams moving at 
different velocities is Characterized by turrulent 
mixing of the streams, accompanied by pressure losses. 
In the course or this mixing an excnange or the momentum 
takes place between the particles UR at different 
velocities, finally resulting in the eguilization of the 
velocity distriputaiors ain the common str Am . The jet 
Ad CH elg er е сосу loses a part of its krrStrTCc ener 
by transmitting it to the slower jet. 


| Ihe loss in total pressure pefore and, after mixing 
ls gua large and positive for the higher-velocit 

jet, an increases with ar increase in tne „amount o 

energy transmitted to the lower , veiocity jet. 
pu Un the resistance coeftricient; тіс неа 
defined as tne ratio of the difference of total pressure 
to the mean dynamic pressure in tae jiven section, wiil 
likewise always be positive., As tə the lower-ve locum 

jet, the eH stored in it increases as 4 result o 

Mixing. The loss in total pressure acd the resistance 
Coefficient can, therefore also nave negative values 
for tie lower-velocity jet [ Ref. 5]: 


Tne projram incorporates this concept at the junction of the 
module ccolinj air and the engine exhaust (if tne system is 


So COALIGULr ed) = The projram assumes the lower velocity jet 
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Figure 2.2 Fan/System Interface. 
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to be the cooliny flcw and the higher velocity jet (о ЕЕ 


exhaust flow. 


Н. EDUCTOR/SSYSTEM INTERFACE 


The eductor discussed in this section is used in the 
engine's exhaust to move cooling air through the сова 
ducting and engine enclosure. ‘there is a mixing of the 
cooling flow and exhaust before it 1S discharged to the 
atmosphere. This section doeS not discuss the eductor 
installed at the exhaust duct exit. The only component of 
interest there is the nozzle as a dynamic loss. The effect 
of the external mixing tube 1S small and can be neglected. 

The module cooling eductor is used on the Oliver Hazard 
Perry caiass frigate. It is shown scnematically in figure 
2.3. The eductor system is illustrated an figure 2. 18 his 
figure snows the geometry and pressure distribution during 
the mixing or primary flow, engine exhaust, and the secon- 
dary flow, поди е сооізпа ELON. А match point concep Cemi 


be develcped for the eductor much like the fan and system 


interface concept shown in figure 2.2. One curve is called 
the gain required and the other tne gain available. These 
curves are shown in figure 2.5. Given the geonetry of the 


mixing area the gain available can pe computed by varying 
the cooling flow while the primary flow, the engine exhaust, 
remains nearly constant ror the desired power setting. Ihe 
gain availatle is a maximum at zero cooling flow. 

The gain required is computed by dividing the system at 
the eductor and is analogous to the system characteristic 


model in figure 2.2. On the downstream side cooling and 
engine exhaust flows move tanrouge the exnaust duct. The 
cooling flow noves through the upstreaa duct. Total PIE 


sure losses can be computed for both aad tne sum is tne gain 


required. Since these computations are taking place at 





Early constant primary flow, engine exhaust, the gain 
uir red at an operating point is a function of the cooling 
flow. The gain required at zero cooling clow is the exhaust 
duct pressure loss under the flow coniition represented by 
the engine exhaust alone. Increasing the cooling flow 
increases the losses іп the exnaust duct and also brings to 
bear losses in the cooling duct. Therefore the required 
gain is a Minimum at zero cooling rziow and increases with 
increasing cooling flcw. 

There must be an intersection oi the gain required curve 
and the gain available curve if the system is to operate. 
шише condition occurs if the gain availabLle at zero cooling 
flow is greater than the gain required at zero cooling flow. 
The intersection must also be far enough totke right to 
provide the minimum cooling reguirement for the load on the 
engine. jhe) Matehing technijue £5 to begin with  scme 
Minimum cooling flow as speciried уу the engine manufacturer 
ӘНІ! March to the right adding a small increment to the 


cooıinğ3 flow until gain required equals gain available. 


I. SYSTEM ANALYSIS 


Sections of the intake ani exhaust ductwork will Бе 
aralyzeaà from node tc node resuiting in the pressure loss 
Bor the section. The sections will be called brancnes. А 
Meme is the starting or ending point of a branch. The 
fittings of a branch will be entered into the program in the 
sequence encountered Бу the flow along a branch. A node is 
an entry, diverging нуе, fan, the Jas turbine en Jine (not to 
Ге confused with the engine enclosure), convergent wye, ог 
an exit. Figure 2.6 snows the six resulting schematic 
representations of a gas turbine installation and the varia- 
wns Of cooling flow available. The numbered dots are the 


nodes. Node 1 is always the main inlet entrance. Node 3 is 
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always the engine. Node ü is always tne cooling far. Node 
6 is always the main exnaust exit. Node 2 may ce Either an 
Eudependentoentry for the cooling flow or the 5ranch loca- 
tion where the cooling flow diverjges from tne combined 
ша е С. Node 5 пау ke either an independent exit Гог tne 
Esolrng rdow or the junction ofr coolinj £iow with the engine 
exhaust. The hashed area is tne engine and tne larger 
rectangle represents the ergine auoduie wnica surrounds the 
Ex: ne anc" 15 a ritting in the €oo0linj flow branch. Ihe 
tranches are designated by the node number at tne begining 
and end of the branch. Ihe reader should refer to the 
user's manual for a complete description of entry of the 
fittings into the progran. 

The system in figure 1.1 would be a class three systen. 
Meenas the cooling flew oranching off the main inlet (11уег- 
gent wye) and joining the main exhaust near the exhaust exit 
plane or the engine (convergent wye). шега со? аза fan 
installec which differentiates it from the class five 


Stem. 
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Ihe Газ1с procedure for system analysis is to assume 
enough Slow and ioss information to proceed with the anal- 
ysis and check the assumptions with continuity of pressure 
at tne nodes with each iteration. iz the pressures do rot 
match, new assumptions are made based on the current 
performance and the iteration is continued until convergence 
is achieved. 

With six different types ог systems to match, six 
differeLt schemes nust be implemented in the computer code 
to handle cverali system matching. Zach scheme must be 
tailored to handie the expected coaponents that make it 
different from any other system. For exanmpie, System six 
has по cooling fan and system one does. System one reeds to 
consider the fan and system interface but system six dces 
not. Appendix A is the complete program listing. Appendix 
E i a flow chart of the most complex system in the 
ГОД гап, system three, and incorporates all possible 


component/system intertaces. 


9. TOTAL PRESSURE GRADIENT 


The total pressure Canan Jes represent tne energy reguire- 
uments of the system. Totai pressure losses in the intake 
and exhaust ducts are inputs to tne engine performance 
Subroutine in the program and are used to determine the 
operating parameters of the engine. Fan and system matching 
is accomplished with the total pressure reguirement. 
Therefore total pressure gradients in the ductwork are most 
important to analysis. Measurement on the other hani 
usually produces tke static pressure Jjradiert. Tres static 


pressure at a point is less than the total pressure at the 


porrt. Figure 2.7 shows a typical representation or the 
pressure changes during riow in a siuple duct. LOSSes in a 
ШІСІ аге due to the irreversible tuuc SDS ibj ое 
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mechanical energy into neat and the losses are used to riot 
the total pressure grade line. Note tnat some rittings such 
as diffusers and contractions cause a Change іт the ӘБ 
pressure guite different from the cuange in total pressure. 
This isa result ora chanje in tne velocity pressure 
tirough a Variable wares, sic eine The sample ргодгащ още ub 
presented in the user's manual, ap ,endix C, саг ое оъеа бо 


produce similar plots cr the pressure Ее ие 
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Figure 





А. GENERAL 


Ihe purpose of the program prepared for this Study is to 
translate the geometry of a gas turbine installation 
Mi luding inlet, exhaust, and cooling ducting inte a one- 
dimensional problem to caiculate tne system's frictional and 
dynamic resistance to air flow and soive tne problem for 
various operatirg conditions. Tne solution will include 
engine performance parameters such as specific fuel consutmp- 
tion, turbine inlet temperature, and mass flow rates. 
Additionally a summary of the duct system performance is 
given by pressure losses for each component and a summary of 
branch losses. Cooling air flow isapredicted py matching 
the system and tne instailed tan or module eductor. 

шррегас че code is utilized for all program “inputs. 
Any number of fittings and combinations of fittings may be 
selected to represent tne user's current design. The system 
in figure 1.1 can be represented by fittings chosen from tne 
menu. About 30 selections from the menu would be required 
to modei the system. The type and number oc selections 
depends on the system's configuration and complexity. Each 
fitting may Lave from one to seven questions posed interac” 
tively to establish the required geometry inputs. With the 
geometry known the program computes areas and coefficients 
necessary to perform the analysis. This data is stored ina 
file called duct data and may be saved for future progran 
runs where ceometry input is not reguired. The oferating 
point is defined upon input of ambient temperature and pres- 
Sure, humidity, horsepower, and power turbine speec. when 
combined with the duct data file the problem may be solved. 


Sis 


В. INTERACTIVE CODE 


Interactive code allows the user to Sit at a comm 
terminal, access a desired program, speciry inputs Еу сура 
at the terminal keyboard, and execute the  grogram. АД 
inputs are requested by statements appearing on the terminal 
screen. Resulting output is written to the user's files 
which may be viewed at the terminal or sent to the printer. 
Tne interactive mode or operation is especially valuable 
because it allows the user, ру modifying selected input 
values, to quickly evaluate the effects of changes to an 
existing or contemplated design. Modification of a system 
is accomplished interactively within tne editor portion of 
the program. The editor offers the ability to change a 
gt to Bg. For example, a mitered rouni elbow could be modi- 
fied to add cascaded turning vanes or a different | elbow 
substituted entirely. Also offered is the avility to add or 
delete a fitting. The addition option does not allow the 
user to add a new first fitting to a branch, however one may 
be added anywhere else. 

The rost lmportant consideration n writing anan r 


tive computer program is what appears on the screen and now 


тар релге; Requests for inputs are in English rather than 
eLjineerirg jargon. Units аге all in the #151151 system 
All lengths are in feet, etc. А41 sogical  cnoices are 


accomplished by entry of one letter, the first ietter of the 
choice. For example, "Y" is tbe reply for yes. АЈ росси 
ChO2Cl 5554355 Gr. A2ndicatec wui ЕЕЕ 
of the question. Shcuid the user not use ore o: tbe с221Сс685 
indicated, the ацеѕііоп wiil Ге repeated uctil ano 
resoorse is given. Default values are avaisaple Гог many 
circumstances to minzmize the Іп Но Ее A default is 
not availiable py simply depressin; the return кет. Тағы 


must elect default values by a озса Е For example 


2 


the Hamilton Standard гі І+ег system installed оп tne 


Beewance Class destroyer 1S avaiiable as a derault for th: 


(D 


сес гісі). Ihe üser SeLicects 61.15 2) апзвегадич al Z ira 
МЕЙІ tc а SUCStTION asking 1£ tne isesi жой] ке to ¿se 


«а 
МЕР Чета rilter system. 


EL OIHER PROGRAM FEATUEES 


ое сэпотаегаяло іп interactive coupuater Eroj-ams 
Ее —.ractice or "user  rroozing'" the inputs. In oraer 
шашы, an zrteractive computer p.ocram shouid rot terainate 
EE utior (1i:.e.,"crasn") if an improper input value is inad- 
verterntly defined by the user. ше пашегтса| валс logical 


iLput two features are incorporated tò protect input tc the 


oO ran. rirst, read statenents are protectel with error 
EC cra of file ietection. А Groblen with irout nere is 
Laniiei ру asking the user to re-enter the value. On 


Mmamerical input if it happens again on tne Same cuestior the 
EXoScram stops execution. s con EE an T Correct numnser 
is properly defined to the program in tne geometry inrut 
рласе, the user is ofzered one last chance to re-enter 
correct fitting data if the user realizes his mistake bezore 
he is asked if he wants to load tne data ror tne fitting. 
ШО Ч5сег 15 assisted here bv a check Cor area cortiruityv 
ОШ оге fitting to the next. Ваша рада іе provided if 
Sent inuitvy is not maintained. Века гос ‘во’ сал и 
Meetirg Erinjs the user cack to tne шепа мат tne rrograrm 
Шу to accept a chcice of fittings Zor usc instead oz the 
erroneously entered fitting. 

The program iS modusarized оу tae extensive use о: 
BeeroOutires. Mose EOM i: tqtescocojgran  I0pt»ve- 
шев Ey allowing tne upjrade and replacement of individual 
ЕО ЦЕ: лез. а р ое 3 ПО СО: 


MESCKS аге usec. Tnerefore сео ааа се үз 


< 


e реет 


ES 


еттш пасы гоп e E ro rana ihe aser nay lecide to change 
the fittings avaiuiazle іг Епе menu, ҒӘР s 1n temm 
code documentation snowS the areas trhat AUST ге с 


accompli Ви 


Atrendix 2 15 а user's па: ча. ñd CITIES The cM 
FIO fal 107 M Э The сагцеі охгізіпз гоя то е а 
t.e rro;ram as installed on tne Naval 2oOSTjlallilte сии 
254 3032 пам. frame ссц мтег айц < 5.3.20 уде поли A 
sampie case 15 lescriped ari sanpile още ме Iovi En š. 
terz2inal session is als гесогісі 52  Á si;c0w tv..cal EM 
Е рыца 


За 


А. GENERAL 


ШООК ЫШЫ ОООО o j о ааш се systexweperlormance се аг 
liar шагіпе ас turbine instailation. аца tie 
Mee рис ШО О х5 roat sucpsections 92 6:3 syster were 
ШЕШ е адс CELL interactior was пеј: естел. m c "ot 
Ee desserrous errors ir tne eSciaatioz of еп ге FESO 
meee buc іс 410 dot rrovide comp.cte  .àicorazation or svsten 
рег- сгоалсе. са СЛЕЙСЕСІСІССІСІ озесооі ив. оч 
ES rot accurate. іше қас Шагсісасагр асассе When tuc 
сеп асі. а госціе есцсгог. 

The rroceas of nanually assicnins a  resistirce coerfi- 
cient to a fitting has been eliminated. МОЯ 1% is possible 
for the computer program to ahaiyze the geometry of most 
ENPtings rapidly and apply tne correct resistance coeffi- 
cients for tne one-dimensional analysis without the user 
looking up any correlations. 

The rrogram rlexibility is denonStrated Бу the ability 
to quickly change input parameters and analyze а system at 
ary operating point. Previous methods anaiyzed comporents 
at full power and then used a proportionality model waere 
losses were proportional to the square of tne engine air 
mass flow rate. This method consistently under-estimates 
duct losses at low power because it does not take into 
account the variaticn or  coolirj <ftiow provided with an 
installed fan or module eductor. At olowwpower the cooling 
flow Can be a Significant contributor to duct losses and the 


previous method can nct predict this contribution. 
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В. LIMITATIONS 


It should be emphasized that any one-dinens foe 
analysis does not hardle riow distortion we Suspected 
rroblems in this area are still kest ijealt with ру (БАМ 
of model studies. Тһе limitation of a one-dimensionai Lodel 
is that a fitting's pressure loss maj ре Known for unifcrm 
flow distribution, Еч is is dizrricalt to ргеацст Ы 


with distorted flow. It is known nowever that tne автор” 


flow situation will have a larger pressure loss, but how 
much is not easily determined. à one-dimensional analysis 
may point to problems with Flow distortion. The program 


recognizes the potential for iow distortion on cerco ii 
fittings such as diffusers and points out this potem m 
Ir a fittirg's pressure 1055 сап vary Significantly with 
distorticn of flow and the one-dia2ensional analysis has 
computed a large pressure loss, the user should flag the 
fitting for futher study by model testing as tke pressure 
loss has probably been underestimated. 

Not all possible duct designs can nave their fittings 
modeled cy the program. Some fittings will be available 
from the program menu and others will be similar to fittings 
listed, rut not exactly. Then there are some which nay not 
be listed at all. I£ the fitting is close, it пау Бейс 
and expectei to give reasonabie results. if tae fitting 
not listed then the user must provide the resistance coefri- 
cient Бу using the "fitting not listed со E Тре data 
for this entry may come from a published correlation GI TTS 
tests performed on Similiar installations. Іс іе іп ББ 
area of correlations where most benefit can ke gained bv 


progran goditication: 
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C. RECOMMENDATIONS 


ШОШ Гг аде ева уе as a stard alore progran, Lut 
some increased utility na; be reaiizei by incorrjoratia, some 
Memthe subroutines in other programs waich would tren izgu*t 
ЕГІ сессронет ад 3PM гесиігелеріс foc ай Operating 
ша Ге instead of pcinrt Ly point user input. 

Тһе General Electric 112200 engine is currently the 
engine within the progran. The engine performance ій tie 
Шо сап is built Бу table interpolation of the puklishel 
performance data. ЕЕ Е Е 155 ое - ‘а or og ran 
Which provides performance uata and it is recoamerd«ed that 
Bus огссгап be sukstituted for tne engine subroutine 
ei rently in the progran. a np eae апу 4сі2%6 
about engine performance rredictions and make the parameters 
more offical. Also the General Electric projraa covers the 
complete perrormance map or the engine wnereas the enjine 
ачольвае. used in this analysis was limitei to 22,507 
horsepower maximum. inere 15 5611. а Little “Sower left 
beyond this value and the program can not currently operate 
Bere. БЕЗЕП шошъгасзе топи сопсегра с the engine is 
iuproving the module temperature out aodel used in the FID? 
Subroutine. The model used produces reasonabie results zut 
is not Lased on test data but on operator experience. 

Пе SU INEPOvenent 1л сгрготгай эесгогаарсе and 
ШЕ су сап ре made by the incorporation of iauaproved fitting 
flow resistance correlations or test lata. 1212105 and 2110 
scale systems should be instrumented to ‚гоу1ае duct pres- 
шие loss data to check the program's analysis. "here tae 
Шо гаш prediction is not accurate new fittinJ7 correlations 
сепа пе developed., Potential fittinjs ror improved пссе15 
are louvers, Silencers, diffusers within iistortei flow, junc- 
mons and wyes (especially where eductor action is desired), 


Eg t oiler tube bundles. йіп Sufrıcient data these 
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fittings could Бе modeled better ава порева 
overall cbjective is to increase bota tne utility anda 


racy of tne program analysis. 
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С. „ВОТЕГ SUBSOUTINE: ІМРПІ5 6ЕОМЕТЕҮ ОЕ 20сТ, OU 201 DUCT Jame C 
Crea eae ТТТЕТТТТТТТЕТТТТІҮТТІТРКТІТТТТІТІТТТІТТІТІТТІЕТ икижкажи аавожиж ак жк анаа С 
С ZO GZT THIS GOING THE DUCT SYSTEM YOU AR? AORKING 4119 ч TO С 
C BE CLASSIFIEO. SYSTEM SUEROUTINE 2OES THIS.  4ITH ТЧЕ Садовае C 
С THE SYSTEM KNOWN, IDENTIFICATION NUMSERS ARE ASSIGNED. THE 1539 С 
C 15 CALLED UP AND FITTINGS ARE SNTZRED зда -HE SISTEM TOA IE с 
C NAMED DUCT DATA. IHZ DUCT DAIA FILE WILL. ЗЕ ЗЕБКІКІЛЛІП 30 с 
С USER WITH A SIX DIGIT NOMBER CF THE USER'S CHOICE. С 
G 
C VARIABLES: WKI AND KR ARZ TRANSPORT ARRAYS ПЗЕР 10 ЕГІ (15 С 
Ç SYSTEM ARRAYS JORKI AND АСАКА. AORKI(NNN,1) ІЗ ТН: C 
C ID NOMBER, AND ЧЯОБКХІ(ЛІН,2) І6 ТНЕ ЕТІТІНС ТІРЕ С 
С JORKA STORES FITTING DATA SUCH AS LENGTHS, AREAS, AND VC 
с RATICS. С 
C 9 * xoxo CXXXI XXXIX О. 
SUBROUTINE 3UILD 
RLAL ЯКа НОНЕЯ 
INTEGER 54525 RKI,ZORKI,7ERM,TYPE,BRANCH,FITIID,GEOCU,DUMMY,4, CLASS 
- DiNENSZON GEOMW (6) “AKIL (2) ака (9) , WORST (20082) . аве аня 
C INST FINDS OUT IF 700 WANT JONG OR 5405 а ж Ванна 
й CALL INST(SORL,TEZRM) 
С SYSTEM CLASSIFIES THE SYSTEM TO ONE OF SIX POSSIBLE SYSTEMS 
CALL SYSTE] 5051, L CLASS) 
C GO TO (lives iss $5) Cu S5 Б 
C GEOM IS THE IDENTIFICATION NUMBES TO BE USED WITH THE FITTING. 
С IZ IS BROKEN UP INTO FOUR PARTS. THE FIRST JIGIT IS THE SYSTEM 
C CLASSIFICATION, 1,2,3,8,5, Un 6. THE NEXT TWO DIGITS ЕН 
С STARTING NODE AND [HE FLUISEING Ч00Е ОЕ ТИНЕ ЗЕЛЕІСІ. ТИЕ ЧЕЗ? 
C DIGIT IS THE FIOM IN THE aRA4CH, 2580 15 COCLING FLOR O DEE 
Ç ENGINE FLOX, T4O0 I15 COMBINED COOLIAG AND ENGINE FLOW. THS IAST 
C THO DiGITS ARE FOR THE ORDER NU4SER OF THE FITTING IN THE BRANCH. 
C EXAMPLE: 113101 SYSTEM ONE, NODS ONE 7C НЕЕ, Звон 
С FIRST FITTING IN BRANCH 
C 
1 ЗЕСИ (1) 5113101 
5564(24-124001 
СЕСИ | 1) 136101 
СЕСЧ (4\ =145001 
PRANCH=4 
CALL FRTCMS ('CLIBRSCRN °) 
ЕСТЬ (0,000) 
CO E | 
2 СЕСЯ (1) 2212201 
СЕСМ (2) :223191 
GECM (3) 2224001 
GECM (a) 2236101 
СЕСМ (5) =245001 
ERANCH=5 
CALL FRTCMS ("CLRSCRN ') 
ARITS (5,600 1) 
GO rO То 
3 37с4 (1312201 
GECM (2)=323101 
GEC (2)=3 24001 
СЕСМ (4) =335101 
GECM (5j = 34 5001 
СССМ (6} =35020 1 
ERANC H=5 
CALL FEICMS ("CLRSCRN *) 
RI; S (G 602) 
GO TO 15 
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акк А АА Аы МАМЫ 


C Evel NG oer OU LINE. ные е Б СНЕ DUCT DATA File (с 
СЕНЕ К Я ЖЖ НИ Жо аки жк яких изя никак С 


С ITE TI. MT Of TIE PROGRAM YOU CAN CHANGE, DELETE WOE ADD à ' € 
E FITTING TO THE TATA FILE. IT 4QULD BE HANDY TO HàVz A COPY OF С 
C ІТ WITH YOU WHEN YOU MAKE THE CHANGES. ALSO THE DATA FILE IS С 
С Е ЕНГ СЕАНСЕ, TO SAVE 4 COBY, MAKE A COPY OF IT UNDER С 
C A CIFEERENT FILE NAME. YOU STILL MUST HAVE A FILE NAMED "DUCT С 
с DATA" TO EDIT. EACH DUCT ВАТА РТЬЕ IS SERIALIZED 3Y “HZ USER C 
C ANC A NEW SERIAL NUMBER CAN 3F ASSIGNED TO THE CHANGED FILE. C 
С THE SERIAL NUMBER APPEARS OF THE COAPUTZD JUTPUT FILE OF SYSTEM C 
С PERFORMANCE. CHANGES ARE MADE BY THE INDEX NUMBER OF THE FIITINGC 
С IN THE DUCT DATA FILE. “HE INDEX NUMBER IS THE NUMBER IN THE С 
C FIRST COLUMN. C 
С С 
С THIS SUBROUTINE DOES NOT CHANGE THE SYSTEM CLASSIFICATION. С 
С TO СЕТ А рІРРЕВЕНТ SYSTEM YOU MUST BUILD IT WITH THE BUILD С 
с PART OF THE PROGRAM. С 
Eben Lice Ue E ou xo ode de oe oe Ro o hE OK oF xA XA XA XA x XX Xx пес 
ЅОРВОЈТІМЕ 5р 
REAL A,4ORKR 
INTEGER Np INDEX, ANS, CHANGE, JELETE,ADD,L,3,5, YES, NC, 3ORKI, P, 2, 
* 
орна нео Я LITE од 
5 2 ; At Z YES И моими 
READ (8,600 SEAL 
DO 10 IS1,N 
READ (8,601) INDEX (I) , NORKI Rn DES &K KI (I, p2) ,WORKR (1a), 
+ &CRKR (I, 2), &ORKR ТОБА 
10 CONTINUE 
REWIND 3 
WRITE (5,592) 
READ | ,903,2N7222, 8 R-22) ANS 
IP ((ANS.EQ.-HANGZ).OR.(ANS.ZQ.JELETZ).OR.(ANS.EQ.ADD)) GO TO 3O 
22 EAIND 5 
WRITE окасы 
GO TO 20 
80 IF(ANS.EQ.CHANGE) GO ТО 40 
ІР ANSE’ DELETE GO 70 80 
Ё IF (ANS.20.ADD) GO TO 150 
C FITTING IS TO BE CHANGED, A NEW FITTING SUBSTITUTED FOR THE OLD 
G WHAT INDEX NUMBER, ü 222 
G 


40 — WRITE. (6,605) 


CAIL READI (M, 5) 
С DO YOU NEED ENU ??? 
50 WRITE (0.996) 

READ | 003, END= 52 ER R= 52) ANS 

I еи TEH: OR. m NS. E EG- NOJ) GO 10 60 

Е v n е е N е б n 
52 т š 

WRITE 15,604) 
GO 79 50 

60 CCNTINUE а 

12 (445. EQ. YES) GO TO 62 

WRITS aff 507 

CALL ВТ (ТТРЕ,5) 

И со 

с 
С CALL THE MENU AND MAKE TEE CHANGE 
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CONTINUE 
GET IHE OPERATING CONDITIONS AND POWER REQUIREMENTS. 
CALL OPZOND(TO,PO,dUMID] | 
ТЯ ЕАН 1S QNGTALZ=D Get FAN CHARACTZR ISTIC 
IF (Class. GT. 4) GO 10 98 
CALL FAN {ROS С, СРМО , СЕММАХ, ОРМАХ ,К) 
CALL PER P j8P 32,20 PO) 
GO To THE SYSTEM SUBROUTINE 231108582 FOR THE SYSTEM 
GO 7O (100,150,200,250,300, 350) , CLASS 
CALL SYs1(SERIÁL,N,4OBRÁI (ORR KY AP, NET, FITIST, 10, PO,HUMID, 
p POE BHOSTD,CFMO ,CFMYAX, DP AAL, K) 
3 
CALL SYS2(SERIAL,N,JORKI,WORKR,HP,NPT,PITIST,TO, РО, HUMIT, 
АРБАП, АСАН, АВИС,АПЯМ 
; SHOSTÉ,Cz34Ó0 ,CFMMAX,DDbMAX,K) 
GO To 400 
CALL SYS3 (SERLIAL,N,#ORKI, FORKR, AP, NPT, FITIST,79,20, HUMID, 
+ ALFAL,ADAB ADHC ADAM, SLEAC АСВ, SCC ACN 
+ RHOSTD, CF40 ,CEMMAX, DPMAX,K} 
со TO 400 
{CALL S¥S4 (SERIAL N,FORKI,4ORKR,HP,NPT,FITIST, 20, 20, HUMID, 
М еле Пао 
5 ВНО5ТА,СУМб,СЕЗДАХ,ОРМАХ,К) 
50 TO 400 
CALL SYSS5(SERIAL,N,WBORKI,WOBKR,HP,NPT,ZITIST,TO,PO,HUAID, 
нос то воо А427 А0"2, 206, А079, ALEAC ACAD, AOAC ACHA) 
{CALL S¥S6 (SERIAL, N, #0 RKI,WORKR, AP, NPT, FITIST,T0,20, HUMID, 
ALFAC,ACHB, ACKC, ACHM) 
"CCNTINUE 
DO YOU WANT TO COMPUTE WITH DIFFERENT OPERATING CONDITIONS 2??? 
ARITZ (6,602) 
READ (5,503, END=820,ERE=820) ANS 
IF ( (383. EQ. Y25). 08. (ANS. 20. 30) ) GO TO 430 
ABITS (6,600) 
A JU eR 
GO T2 195 
CCNTINÜE 
ІР (А85-20.153) GO 72 95 
FOEMAI(IG,/, 13 
FORMAT 13,3k,16,3X,12,3X,7 10-4, 31,7 F10.4 3X, F 10. 0, 3X ,P 10, 0) 
FORMAT(* “DO YOU WANT TO cOMPUTE WITH DIFFERENT OPERATING CONDITIO 
5 (YZ EN 
FORMAT 1) 
FORMAT ' YOU MUST ENTER A LETTER INDICATED IN THE BRACRETS."') 
«ае А 1 
END 
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ВТО SSL CRED TEN SEND 


Cx xod UR OUO XA X KOR OK KOR R жижжиж Xo жж жж жж яя ниж жк ин жи жк C 


С FIIs SUBROUTINE CALI S LOAD à SUBSOUIINE TUAT TBANSPERS ТН? С 
г Ша A А ща Aa Rc usb = 
С ATA OF A PITTING TO THE SYSTEM STORAGE ARRAYS 4OPKI AND ZORKR C 
C ex Rua aRAREREARA RXRAXREAARER SER ARERRAXAS RRSRRRSEE AF GXKG RES ERRAT KK HK Kk 2 * C 
SUBEO UTINE S SELECT (М, SORL,GEOM, TYPE, KORKI, NORKR) 
4 А 
INTEGER Ў акт, WORKI,SCRL,GEOM,TYPE 
DIMENSION зоа КІ (200,2), AORKR (450 4), SKI (2),3 „ЯКЕ 3) 
с CHANGING THE NUMBER OF FITTINGS BEQUIRE TC THE FOLLOWING 
© СО ДО STATEMENT AND THE ADDITION Or Aã CALL ти THE SUBROUTINE 
ЩЕ pe А E ل‎ ° 
CERE 0o 576 748.9. 1911. 12-13, 14,15, 16, 17, 18, 197 20 
+ 954 ань 252527, 28 уә 2222 Te? 
1 CALL 71701 (5084 GEOM AKÍ, RKA) 
CAIL LOAD (A, CECU, IKI, AKR, HORKI, HORKR) 
4 
2 CALL FITO2 |5051 GEOM, EKI, AKR) 
CAIL LOAD (M,GEOM,WKI ,4KR, 4ORKI,WORKRE) 
бо То 40 
3 CADDSEPPOS (59 ЕТ GEOM , RKI,RKR) 
CAIL LORD (N, SECS, aKL, ¢KR,WORKI, AOBKR) 
со то ч 
4 CALL ЕТТОЧ [SORI GEOM, HKI, WER) 
CALL LOAD (M,GEOM, AKI , WKB, WORKI, WORKR) , 
5 CAIL FITOS E EO Я, ЯКІ, Яка) | 
Cail LOAD (À, GECA, HKI, WKB, WORKI, WORKA) 
Š CAIL E SORL,GEOM,WKI,WKR) 
АРОК (M-GECN, AKL . KB, NÜRKL, ORKR) 
GO TO 40 
7 C ALL ў. SORL,GZOM,WKi, WKR 
саги позор. (M, GPOMPARIOSKHANORKI,4OBER) 
GO TO 40 
8 CALL FITO8 (50БІ,СЕОМ,ЯКІ,ЯКВ) 
CALL 1040 (Ч, СЕСЧ, ЯКГ ‚ ЯКЕ, НОВКТ , ЧОВКН) 
9 Gall FITO9 (SORI,GZCM,WKI,RWKR 
CALL LOAD (M,SECM,a4KI,BNKB,WORKI,WORKR) 
GO TO 40 
ШІ СІШ ЕІ:10 (SORL,GEOM, (EI, E) 
CAII. LOAD ( , СЕСЧ, ЯКІ, КВ, ЯОВКІ, 4ОВКЕ) 
11 C I (SORI GEOM, RKI, (яка) 
САЦ 10 AD ( СЕМ ZEPOWKB. WORRI,AORKB) 
TO 4 
12  CAIL FII12 (SORL,GEOM,WKI,WKR) 
GAIL LOAD (M,GEOM, WKI , ¥KR, WOKKI, AORKR) 
GO To 40 
ШІ | СКІ TITY (SORI стон, WKI, WKR) 
CAIL LOA MM GEC ETI UKE KORKI WORKER) 
we «GATE FITi4 (SORL,GEOM, AKT, HR) 
CAIL OAD ( SG EC S AKI , RKR, WORKI, HORKB) 
ILES TT 5S (SORL,GEOM, “ЯКІ, Яка), 
CALL 2040 ( CECH a АКТ ЯКЕ НОБКТ,НОРК В) 
16 CALL FIT16 (5051 СЕОМ,ИКТ, ЯКВ) 
CALL ІСАр (М4,6Е00,ЯКІ,ЯКВ,ЯОЗКІ,ЯОНКВЕ) 
GO TO 40 
17 CALL FIT17 (SORL,GEOM,WKI, WK а). 
CAIL LORD (4, GECM, «КІ, ИКЕ, ЯОБКІ,ЯОВКЕ) 
19° CAIL FII18 (SORL GEOM, “ЯКІ, Яка) 
CALL LOAD ($,GEOM,WKI,WKR,WORKI,W4ORKR) 
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ETITITDIJGEUDQCELDEONAMOSOOTHCRBADIUS, ROUND CROSS-SECTICN C 
Ска ЖЖ ЖЕНЯ К К НЯ И EE EEL EE RETR ERE EERE RE EEC 


Е REF. ASHRAZ HANDBOOK, PAGZ 33.33, TABLE B-3, FITTING 3-1 C 
C CURVE FIT 70 THE TABULATED DATA C 
C SHCRT FITTING, FRICTION LOSSES NOT INCLUDED, CCNNECTING DUCTS С 
С TO THE CENTER CF IHIS FITTING. с 
ЕЕ ik 

SUBROUTINE PIT03(SCRL,GECM, “КІ, ЯК В 

о а ИДЕ 

DIMENSION ЯКЕ(2),4КЕ (8) 

45175 (6 600) 

CAIL 52408 | ,5) 

58175 (6 601) 

CAIL BADR | 25) 

RRITE (6, 60 ) 

САК. ЕЭВЕ ы» 

=U. < 

KI32T7T3:20.0306 *18ETA**0,7825 

CPRIME-0.029469*EYP (2. 5627*(1.57138-(R/D))) «0. 11746 

С КТНЕТА “С РЯТМЕ 

ЯКТ(1) =СЕОМ 

ЫКТ (2} =3 

ЯКЕ (1) =АВЕА 

WKR(2)20.0 

НЕТ 
600 FORMAT('. YOU HAVE SELECTED A SMOOTE RADIUS &RÓUND CROSS-SECTION'/ 

+ . М у = 5 < ша 

: PEOR ТИ **FIRST QUESTION, WHAT IS THE CROSS-SECTICN DIAMETER 

v 2 
601 FORMATU ЕНЕВ ОТНЕ BADIUS OP THE TURN OP THE ZLBCZ MEASURED TO* / 

+) (THE CENTERLINE CF THE DOCT.*) 
602 FORMAT(’ LAST CUESTION, ENTER THE ANGLE OF THE ELBOW TURN. (DEGA 

EES 1 

* = 

RETURN 
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Cx ox иж CURRO OX КВН В а и жи как жж С 
с EF. ASHE AE I ANDBOOK PACE 33.33, TABLE 9-3. FITT e e Č 
с CURVES FIT THE TABULATED DATA FOR SACH NOMBEE OF SEGMENTS. C 
с THIS IS A SBORT “PITTING, FRICTION LOSSES NOT INCLUDED, MEASURE С 
С CCNNECTING DUCTS TO THE CENTER Or THIS FIITING. С 
Qo x XA AX X Xx kt Tk x XO жжж ххх хуи хх KEKE KX жижиж жи иж #ж жи жжж и «жи »« Kk OC 

SUSROUTINE ЕТТОЧ (5081, = 
REAL D IR, IRB 
INTEGER 5081 СЕСМ ЯКТ, Н 
DIMENSION HWKÍ (2), KR (0) 
5 ЯВТТЕ (6 ӨЗІ 
САТІ, БЕКОВ (4,5) 
zz 443-2223); 6R- (5. 87. 5) GO TO 5 
WR == (6.601) 
САГ 254281 ,5) 
RRITE (6 60 ) 
GAL БАРЕ | 5) 
АВЕЛ=О. 78 54*6 **2 
M2N-2 
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GO TO 
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GO TO H 
30 С= 1. 0456 жЕХРр (1.743 13* (0.01219-8 2) ) +0. 15776 
40 CONTINUE 
ЯКТ (1) =ЗЕОМ 
ИКТ (2) =4 
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YO ete MEASUREÉ ТО ТИНЕ CENTERLINE OF THE DUCT?*) 
5210 N 
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CALL TABLE (T1,X, XOUT , CPRIME) 

CALL TABLE (T2,X,XOUT,C2RIME) 

GO TO 59 

CALL TABLE (T3,X,XOUT,CPRIME) 

CONTINUE 

TE 100: (0). eT. ОО вас я 2ў аза аўў GO 70 60 

WRITE 

GO TO ($5 

C=CPRIME*KTHETA 

SKI(1)2GEOM 

ИКТ 12) <3 

ЯКВ (1) =АВЕА 

ИКЕ (2) =0.0 

ЯКР (31 =С 

ИКЕ (Ц) =АВЕ 

FORMAT(* YOU HAVE SELECTED A SMOOTH RADIUS RECTANGULAR ELBOW WITH 
+ SELIZTERS.'/' IT ҸАҮ HAVE 1, 2, On 3 SPLITTERS. ') 

FORMAT (|| *#FIRST QUESTION, HÓW MANY SPLITTERS ARE IS THE ELBOW (1 
252 : 

Ат і “HAT IS THE HEIGHT OF THE EL 30H?" / 
+! (THE CROSS-SECTIONAL DIMENSION PARALLEL TC THE TURN AXIS)'‘) 
FORHAT(' ИНАТ IS THE WIDTH OF THE ELBOW (THE CROSS-SECTIONAL" / 
e! DIMENSION IN THE PLANE OF THE 2033) 2*] 

FORMAT(" ЯНАТ 15 ТНЕ RADIUS OF THE SLBOW, MEASURED TO THE CZNTER' 
+/' OF THE ELEOW CROSS-SECTION?) 

FORMAT( LAST. QUESTION, WHAT 15 THE ANGLE OF THE TURN (0-90 DEGR 
+EES) 23 

ORM ATY® CROSS-SECTION EXTREMELY NARROW, RE-ENTER BETTER DATA.') 
RETURN 
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ІТТІМС 11: ТІРОН5 90 JEGAZE RECTANGULAR IN 2e p R n C 
Сина авовивосвоакосково Вазсиаванааа вохизт ана с иззеванеснасввк аежиж а 
C НЕГ. ASHRAE HANDBOOK, PASE 33.32, ЗАВГЕ В-3, О О C 
C СУЕТЕ FIT TO TEE TABLE cATA C 
с C 
Сив авнаннанавекавиянена ваен вва вананаавнаненнанеен ка ватаякаян ания какакаС 

SUSROUT INE PITTI soa [gto АКЕ ОКК) 
REAL WKR,C,AREA ЕБІМЕ Х,Ү,К 
INTEGER 56 L Е ркы 
DIMENSION В 
WRITE (6,600) 
CALE 32408 (4,5 
WRITE (6,60 Me 
GALL BEADR (8,5 
WRITE (8, 60 " 
CALL SEADR ( 
X=L/H 
TF AX GT ; AND (X 22877 60110 10 
Á“ + Ф е е е Т е ) 5 
о : 85 
GO TO 29 
10 c^ 1 49- Е „21052) #х+9. 1399) #х-2.168) #Х+0. 0545 
20 | Co! 
K=O. 4 704 85 x2 (- -0.3558* 1) «0.67 
СРЕТМЕ<С #К 
2822.9» [H* 4) И (Ней) 
ААЕА=Н® 
ЫКТ (1) =СЕОМ 
Е 
ЯКЗ (11 =АВЕА 
ЯКВ (2) =ОН 
ЯКЗ (3) =ОРЕТЧЕ 
AK RIGVSAREA 
600 FORMAT(' YOU HAVE SELECTED A SERIES 90 DEGREE RECTANGULAR EL304'/ 
+ SET IN A Z-SHAPZD CONFIGUAATION.'/* жертват QUSSTION wna Tu 
tHE HEIST ӘР tHe БЕВОН CAOSS=S=CTEON? I (DIMENSION IN THE PLA 
+ НЕ THE T 
60 1 FORMAT QHA T IS THE WIDTH OF THE ELBO3 CROSS-SECTION?'/ 
+! (2 E DIMENSION PARALLEL TO THE TURN ARIS"). 
602  FORMÁT(* LAST CUZSTION, “HAT IS THS LENGTH BETAZEN CENTERLINES*/ 
+' OF THE "Z" ENTRANCE AND "а" &411?') 
RETURN 
END 
T3 
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C REF. ASHRAE HANDBOOK, PAGE 33.33, TABLE 8-3, FITTING 3-12 = 
C CURVE FIT TO THE ТАЗОТАТЕО SATA C 
C C 
ИИ aaa 


SIERO ULINE fl. iz (SORL,GEON, AKL, NAR) 


АЕ С AREA DS bea, od, CPRIME,X,Y, K 
LILE ER SORL SECM, IKI 
DIMENSION ЯКЕ (2) яка (4) 
&BITE (6,600) 
CALL SEADR (H, 5) 
AR ITE (6, 00 1) 
CALL EADR | ,5) 
WRITE (6, 002) 
CALL aEADR (L,5) 
X=L/4 
У=ЯИН 
TEE 0.0) 20. (1.17. 1.4)) 60 то 10 
EDs. | sa AND ste LT. 2.0 | GU T9 20 
TEU B2 O) ANDOMNIGII.N.0)) GO TO 30 
Cz «40. 10 *X 
CO TORVI 
10 C^ ((((15 79380391) 75.07 366) X«3.5357) «1*2. 29886) *X« 1. 20 
(5 
20 1665. 57198-50118) +5 1.0057 | "mi 
3 e - 
30 C7 (140.30983*X)-0.206 799) *X«1. 154425) *X «1. 702965 
O — CCNTINUE 
K20.4 7/08 *E XP (-0.3558* Y) «0.67 
СРЕТМЕ=С *К 
5822; je [H* V) И (нея) 
AP EA=H* 
WKI(1)=GEOM 
WKI (2) =12 
ЯКЕ(1)-АЗЕА 
«КЕ (21 =ОН 
ЫКЕ (3) =СРВТМЕ 
f ЯКА (4) зАКЕА 
600 „БОВМАТ( 100 HAVE SZLECTED A SzT CF 90 DEGREE RECTANGULAR £13045 
+IN DIFFERENT PLANES.'/' **FIRST QUESTION, WHAT IS THE HEIGHT OF T 
tHE ELEQW CROSS~SECTION?'/* (2IXENSION ^N THE PLANE OF THE TURN 
+ 
601 FORMAT ( qHàT IS THE WIDTH ОҒ THE SLBOW ZROSS-SECTION?'/ 
+ (THE DIMENSION PARALLEL 20 IHE TURN ALISIN) 
602 FORMAT(" LAST QUESTION, UHAT IS TEE LENGTH BETJEEN CENTERLINES'/ 
+‘ OF THE "Z" ENTRANCE AND "Z" EXIT?') 
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C FITTING 16: MAIN SECTICN CONVZREGING ЧЕ 

жк ж ажа ал жала 59 На 2 зв икон ижи ана зо нни о ево ножи е ниянизС 
о HANDBOOK OF HYDRAULIC RESISTANCE, SECTION SEVEN 
PAGE - 


AAA 
ССС 


C*****5*ts5*X5*5**4&* *ikckck co x A 8 0 ааа A 8 G0 G8 фах та ЖҰЖ ааа * C 


о FIT16(SORL,SZCM,WKI, AKR) 
НЕДІ АКВ, АМ 


THTEGER SORL £ GEC, ЯКІ 

DIMENSION QKI (2), KK (4) 

8a ITE (6,600) 

CALL READR (ам, 5) 

ЯКІ(1)56Е0: 

WKI(2)=16 

НКЗ (1) =А М 

WKR(2) 20.0 

яка (3) 20. 0 

ЯКЕ (4 | зам 

600 РОВМАТ(' YOU HAVE SELECTED THE MAIN SECTION OF A CCNVEBGING'/ 

+ ЧҮ, THE ENGINE EXHAUST ALON® SHOULD 3E FLCWING THROUGH '/ 
** THIS SECTION. IT SHOULD B82 THE LAST FITTING CF THE cBANCH.'/ 
e't  **JUST ONE QUESTION, ЯНАТ IS THE CROSS-SECTIONAL AREA OF THE'/ 
+ MAIN ERANCH?*) 

РЕТПЕН 

END 
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C EF. IDEL'CHIK, HAJDBOOK OF dY23AULIC &ZSISTANCE, SEGDRONMEUIN NER: 
С e C 
C С 
СОР ыма е ха как Жк ктк кек Ез «жж seuss 
SJEROUTINE PITIS (SORL SECH, RKI, AKR) 
3 iL GRS 30 Br ,K1,K2, A0,a1,ALTA,BETA, THETA, CELE CE R EW 
INTEGER 56580 SoRL,8KI,ÁNS,TES,NO 
JIMENSION WRI (2), ARR 44 
JATA (55/4 187, МОИ" У! 
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C REF. ASTRAE HANDBOCK, PAGZ 33.22, TABLE B-7, FIZTING 7-8 Č 
С CURVE FIT TO TAEJLATED SATA, BASSO ON O0CS SPER AND PEN с 
C FREE FLOW AREA. С 
255556 ЕЕЕ жи я аго тата ЕХ ЖЕЖ кк жж «жж ЖЕ s.s... a... И 
SJPROUTINE РІТ23 (SOR Lr GEO, AKI, KE) 
REAL „КЕ, РОСТА SCRNA,N 
ТІТЕСЕР 508 ИС сечи 
DIMENSION sWKÍ(2),4KR(4) 
481756 500) 
CALL BADE | UCTA,S) 
WRITE (6 60 | 
CALL 24254 C& NÀ, 5) 
Y= SCR NAZDU а Б 
T KC 902 1+н) -92. 445) #+32.066) #Н- 1. 9557) #н+0.025 
ЯК =С РОМ 
KID 
НКЯ (1) DUCTA 
КЕ (2) 20.90 
яка (3) +С 
“Ка (4) <ВИСТА 
600 FORMAT(' TOU HAVE SELECTED A SCREEN OBSTRUCTION IN THE DUCT.'/ 
+ ' а ЖЕТВ5Т CUESTION, AKAI IS THE DUC C05 SECT T | 
601 РОВМАТ(' LAST CUESTION, WHAT TS THE FAZE FLOW AREA OF IHE SCREEN? 
+ ç : 
RE TUR 
END 
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ЕТО Чо 24 ГО OVER SSN TRANCE 


С С 
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C 
C CURVE FIT TO DYNAMIC LOSS IJFORMATION, 
S 
ОСЫ r sss ss 

SUEROUTINE FIT 

REAL Dr, 0B,N, 30C 

INTEGER AKI, SORI. GLO 

CI*ENSIZN ЯКІ (2), AKR (4) 

WRITE 32041 

САТЕ READE (DX, 5) 

ий мес (6,6 1) 

CALL 52355 0 MED 

WEITE (5,5 2) 

CALL a£ADR (4,5) 

SRITZ {е 993). 

CALL aEADR(DUCTA,5) 

F=N®DB/DX 

С-62.143чЕХР(-4.27543%чЕ) 

4KI(l)sGEOM 

WK I (2) =2 4 

Fl) DUCTA 

WKE(2)20.0 

ЯКЕ (Зў С 

ЯКВ (3) =БИСТА 
600 FORMAT(' YOU HAVE SELEC 

+ t $*FISST QUESTI 

+ ' LOOVER OPENING 
601 FOPMAT(' WHAT IS THE DI 
| + От DISTAR 
602 БОЯНА, JJW SANY CPENI 
603  FOBRSAT(' LASI ;UESIION 

* ' INSIDE THE Log 

RETURN 

END 
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C COMPOSITES LOSS COEFFICIENT SASZD ONCASSUUDENSMGONDNONGOISP 
Fa4LEION AND А 5000ЕЧ EXPANSION 
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«ак TK TK kk kk k; kk акк их EAE ж 


50ЕЗ097 INE Е NS (ЯКЕ) 


m 


morb «Ка 5 H CX a Al, DR, АСТ СЕД СЗ ИЕН 
INTEGER SOL S tC 
Aa aL G 
CIMZNSION ЕНЕ а 
aRIZE (0,600) 
САИ FEADR (5,5) 
WRITE (6 o0 1) 
CALL 35408 (2,5) 
GRITE (6,00 ) 
C Dn 354051 ,5) 
&RITE (6,603) 
CALs mHEAQR (d, 5) 
SRIIE(b 934) 
CALL 12008 £ 1,5) 
4RITZ (6,00 } 
CALL ЕД ОВ (Ч,5) 
А0=СХ*Н 
А1=Ч*С*Н | 
. 0Hz2, 0*3 *H/ (G* B) 
с ы ен 
шш ` + E 
C120. 114 1* ((R/DHe0. 1) **14,0425) *(1.0-A1/42) 
C FRICTION 
С ба ада ШІ са 
С3:0.274|1,2-4%07с1)%е2%0.02 
C CCMPOSIZT2 COEFICIENT 
C=C1#C2+C3 
SKI(1)23EOM 
#КТ (2) <26 
WKE(1)=A0 
ЯКЯ (2) =9.0 
RAE 3 i 
AKA (U) =A 
600 FCBMAT(' YOU HAVE SELECTED i MULTI-3AFFLE TYPE SILENCE 
+ ' | EACH BAFFLE HAS A STREAMLINED SHADE. Е 
+ ' USSD IN HE INLETS OF TH2 DD963.'/ 
і + t **PIBRST QUESTION, ЯНАТ IS ТНЕ САР ВЕТИБЕМ ТНЕ 
621 FORMAT(' ЯНАТ IS “НЕ THICKNESS OF THE 3AFFLZS?') 
602 FOEMAT(* WHAT IS THE LENGTH OP THE E 
603  FORMAT(' WHAT IS THE DIMENSION OF THE 3AFFLES 2AFALLE! 
э 
к; 
o 25 4 e АА iiu. ‘ as i . 2 л) ب‎ 
50% FORMAT (* WHAT IS THE DIMENSION OP THE MAIN DUCT ACAOSS 
+ 
0492 Ба 2 1 N, А 5 Ë e ae 
595 ropar ' LAST QUESTICN, HO4 MANY GAPS АВЕ THERE?!) 
awe Н 
END 
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C prottug 27: CA URBINS MOLILE C 
Стжжжжажжжж ох ких Хи це иж жи жк изжажа и ижи а жука а 
Е REF. SENZBAL ELECTRIC САТА „105525 IN THE MODULE S3ASED ON THE С 
E MASS FLOW THROUGH 212 OLE VI CITING ASD HER. IS C 
C SUBROUTINE JUST . 3:2 THE MODULE. NOT THAT “HE *ODULE IS C 
С THE FLOm PATH АСТ Баз ÉNGINZ. osos ILL 3 ІН Та: СООТІНС РТО. 
CeERRKKKESEE EEE ES ее. жж) жж її їнїнїн» С 

SUBROUTINE FIT27(SORL,SEOM,8KI, AKS) 

RON АЮ 

ТУПЕЕ Са БОР ЗЕЕ RT 

DIMENSION як! (2) „Яка (4) 

ARITE (6,690) 

ЯКТ (1) +ЗЕОЧ 

КТ (21-07 

ЯКВ (1) =1.9 

ЯКЕ | 2) 1.9 

ЯКЕ (3) =1.0 

икн | 4) 51:09 
600 FORMAT('* YOU HAVE SELECTED THE GAS TURBINE MODULE AS À PART OF'/ 

© ТИЕ GCOLING FLOM PASSAG?Z. NO QUESTIONS, JUST NEPDED'/ 
+ ' TO KNOW WHERE YOU WANTED THE м9091Е."“) 
RETURN 
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С DUCI DATA FILE CUT2UT SUSROUTINE C 
C C 
C ЯБІТЕЗ THE SYSTEM ARRAYS XORKI AND SORKR TD THE DUCT DATA FILS. С 
С ALLOWS THE USERITO SERIALIZE БАСА EI E CRECEN Č 
C S4AFNINS* ARITES OVER OLD FILES, SAVE THEM JUNDER A LIFEBENT NASZ.C 
Сажин аки ааа S32 S222552*239955295 S455995592329924 3925 S222 а С 
SUBROUTINE SU1CUT (NORKI, 3ORKa, M) 
REAL 4OBKA 
INTEGER WORKI,#&,SERIAL 
DIMENSION RORKI(200,2) ,NORKR(200,u) 
ARITE (8,600) 
GATES 15215 RIAL,S) 
ЯВТТЕ (3,6 1) SERÍA 
WRITE (8,602) М 
ВО ES 21175 603) I, ZORKI (I, 1) ,RORKZ (1,2) ,RORKR (Lp UOI E 
IRITE е ) a І (1 я - ), e О ), 
+ WORK RE (1,3), sOKAR (I, 4) 
10 CONTINUE 
REWIND 8 
60) PORMAT(’ WHAT SERIAL NUMBER WOULD YOU LIKE TO GIVE THIS DUCT DATA 
+ ТІГ22%/б You MAY USZ UP TO A SIX DIGI? INTEGER NUMBEB. ") 
601 FORMA? (I6 
602 БОБМА? (13 
603 EOENAT 13,3X,16,31,:2,3X,F 10.3,3X,F7 10.42 ЭХ, РО, nr О 
? - 
СИПЕТЕ захж аа ава « «за «# «в # ### # «# # # # «а ## « k C 
C RIAL NUMB3IR AZAL SUSROUTINT; PRSS FORMAT C 
C PREVENTS THE INADVERTENT ENTRY OP NOLL DAT {HITTING THE RETURN С 
C KEY 8124 NO ENTSI) AND INCORRECT DATA, 2413 0002582152055. Е 
Е TT ALLOWS FREZ FORMAT IWEGT. 240 NULLS XILLS таз 2ЕОСНАЙ: 
(***5****S5*X*X5€"5*x€"XX5*59*5*555**5*5* ЖЕКЕ ЖЕ ЖҰ ЖҰЖ ЖЕ ЕЕ ЕКЕ жж Ж 552552255240 
SUEROUTZNE REALS (ANSR,FD) 
REAL ANSR 
INIEG:R COUNT, EL 
COUNT =0 
10 CCNTINUZ 
COUNT=COUN T+ 
Ir (CCUNT.£7.3) GO [0 20 
CALL FRICMS (*CLRSCRYN °) 
aRITE (6,600) 
бй то 09 
22 CONS NUDE | 
Read (FD, *,£ND=30,2R8=30) ANSE 
БІТІР? 
30 BERINL ЕО 
WRITE (6, 001) 
зе те то 
40 CONÎ INTE 
600 FORMAT (////* PROGRAM RILLED - T30 NULL STRINGS ENTERZI!'/) 
601 204420 WARNING: NULL STRINGS ARE NOT ALLOWED, ENT2S UMESIC 
«АТ VALOE. 
END 


Loe 


C******5555 555 X52 аквожан ии жи кк жк ии аа кик к каки аажкии EERE eee C 
С INTEGER READ SUZROULINE: FEES FORMAT _ с 
С PREVENTS THE INADVERTENT ЗЧТЗУ OF NULL DATA JEITTING SET C 
C KEY atta NO ENTRY) AND INCORRECT DATA, THIS КСОТІМЕ I5 USED. С 
C RE ess een анна та а аа чу 517216142 2505519. С 
ӘТЛІЗІТТІТТТТГЕТЕЗІІЗ ІТ 522325233223 2225225425 522322222224 22227%777% 
SIE OUI SALI (4S 7.72) 
ІМТЕСЕЗ -ІПЧТ, АМОІ,?П 
CO NTD 
10 CONTINUE 
CCUNT=COUNT+1 
ТЕ CONT TI GO TO 20 
CALL FRTCMS (*CLRSCRN ') 
WRITE (6, 600) 
GO TO чо 
20 CCNTINUE 
RZAD (FB, * 252 230, E88 230) ANSI 
FETOR 
30 REWIND FD 
ARITE(6,60 1) 
GO 70 10 
40 CENT INUE 
600 FORMAT (////' _EROGRAM XILLED - TWO NULL STRINGS ENTZEED! ү! 
691 FORMAT S ots аа А бае аба ARG NOT ALLOWED, ENTER UMERIC 
EAS UIS 
EN "E 
НЫ ках ких E E E RE f Jk Û f Û жС 
В OPERATING CONDITIONS INPUT SUBROUTISE: TEMP, PRESS, HUMIDITY C 
С SUBROUTINE TO INPUT THE FOLLOWINGS DAEAMETERS: С 
С AMBIZNT TEMPERATURE (2552245 ? C 
С AMBIENT PRESSURE (PS à Ç 
c RELATIVE SUXIDICY (GRAINS PZB POUND 4ASS DRY AI2) C 
C 
ЕЕЕ ЕЕЕ СІ: 25 ЕО 2G) A@Gs2T SELATIVE HUMISITY (%) AND C 
са сара аса аа 58020152 ТАТСЕ ІҢ ТҢЕ 2506614. THAT С 
Ge MCDISICATION SHOULS SE ACCCYPLISH=D 4296. C 
C***t**s55*t* *oxod «а ### ## ##а E #а# # «а##### жж ЕЖ ЕЗ ЖЕ «ЕЖ ЖЖЖ ЕЖ жж 555552252543 546 
SUERCITINE OPCCND(TO, 20, HU HID) 
RETO 20 ОМІ 
ARITE (6 909) 
CALL 22408 { 0,5) 
ЯЕТТЕ (6,60 ) 
CAIL 354944 0,5) 
WRITE (0, 60 2) 
CALL HEADS (HUMID, 5) 
АГ ЕСН (2 CLAS ÇEN 
600 FORMAT(' THIS PORTION OF THE PROGRAM INPUTS THE ENVIRCMENTAL COND 
| -ITICNS.'/'  d4BAT IS THE iMBIEZNT :ZMPIFATCUSE (2EGREES ©) 2°) 
601 FORRAS I, WHAT iS ТЧР AMBIENT PRESSUEE (PSIA)?"') 
592 228722 ШЕННІТ ІС чо СЕТТЕ 117450ІТҮ (254125 CER FOUND AIR)?') 
ENÊ 
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фах &$593* хх Еж ж ж изо киа а с 
POTER 22INT INPUT SUB BCUZINE (HCaSERCHEE, PCHER TUPBINE СРЕЗО! С 
DCES A PRELIMINARY TEST TO INSURZ 75520: LIMITS APE WOT EXEC БОШ 
POR THE JZERATING PCINT 850955155. 955% І5 420025250 TO єс сан 
A CIFTERENT OPESATING POINT го зат TIS НЕ N C 
за а ха а «я ея я какая ии каки r жС 
SUBROUTINE 288 PT(H2 327 1) 221 
AZAL На, ЧРЕ Ро, 20 гасо раса весла ае NA 1565 | 
48122 (6 555) 
CAIL AZA DR (H , €) 
481ТЕ (6,6 1) 
CALL ЭВА 58 (МРТ s) i 
Ir ((NPT.SE-.1290.0) .AND. (NPT. 22. 369020) REO TOME 

«RIT? (6,603) 

225107720 
CONTINUE А й 
iptg2(5177.5-7.0527*10)*(8. 3275-0.0012* 70) * N27 
WRITSIO.002) 
GJ TC 10 
CCNTINUZ 
OO. PRICSS('CLPSCEN зи 
FORMAT(' INPUT THE 2OWER SETTING YOU DESIRE.'/ 

© sad IAT IS THE jORSPPOHE RON 
ТОБМАТ(' **4HAT IS THE POWza TURBINE SPEED (£23) 2") i 
FORMAT’? HORSEPOWER IS NOT ON IdE 2EZEFORMANCZ ИАР, PICK 3 LOWER 5 
CRSEPC ER.) ш 
FORMAT(' OWNER TURSINE RPM LS NOT REASCNABLZ. IT SHCULO BE'/ 
' 1200 10 3600 REM.  &RE-ZNTER.') 

RETURN 
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COMPUTE OUTPUT SUBROUTINE: PEINTS SYSTEM DATA С 
THIS SUBROUTINE WRITES TO THE OUTPUT FILZ. IF YOU HAVE AN OUTPUTC 
FILE ALREADY ІТ WILL BE WRIITEN OVER BY "ЧТЗ PRCGRA*. I= YOU C 
NANT ZO SAME TEE PREVIOUS 4250175, АЕЧАЯЕ TES FILZ. IF YOU 400 С 
OR CHANG? FITTINGS YOU MUST 4AKZ SO4E CHANGES HERE. С 

* ooo G0 X 0A Xx 0 AG AE RO oo XC GR G0 GR X COR XX A 0 XA X X 0X XGA JG GA xA xA X G x G X A GG 0x Xx Xx xA Xx **k*km*5x****5*X*x(c 
SUBROUTINE JUTPUTI(Z0,PO,HUMIO,HP,NPT,N,SOEKI,DP,FITPY,INLOSS, 

+ EXLCSS, IC, H2, 43,29,79, SFC LS NOI SERIAL 
REAL 10,20 ,НОМТЕ, Че, МЕТ, ОР.РТ-В\, -16055, 25.059, ИЕ ры 
+ Я8,Р8, 18, 520, 1504, 16 
INTEGER N,sORKI,SERIAL,T (PEZ 
DIMENSION 28 (209 FITPY(200) BO RKT (200, 2) 
WeITE (4,600) SERIAL, TOPO MM ID IB BE 
ЯВІТС (37601) ІЗІ,055; ЕХІ055 
WRITE (4,602) WC,42,4 8,28,T8,SFC,T54, NG,TMOD 
20 100 121,5 
АТЫН УК МАН 
fa 2, ls 
Біты (ПЕНИ Из, 19636 
со T9 (1, 2239 25 ТО То ая 15, 16,17,18; 1922/0) 
а ааа 541252627 29^ 287 2$ 130] 
&RIZZ (4,503) W НОРКИ AOKI (û seti. ЕТТРУ (1) 
aRITE (4,604) 3ORKI(I,I) -WORKI (I, 2) ,9P(I) -FITPV(T) 
3 - 
Я8115 |44605) RORKI(1,1),WORKI (1,2), 0R(L ШЕШІ 
4 
WRITE (4,606) WORKI (1,1) ,3ORKI (1,2) , DP (1) ,FIT?Y(I) 
GO TO 100 
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CONT 
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GO 20 
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GO TO 100 Е 
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GO TO 106 
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GOT 
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GO Т 
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а E 
ARIDE (4,620) WORKI(I,1),WORKI (1,2) , DP(I) ,FITPV(2) 
RITE (4,621) KORKI (1,1) , NORKI (I,2), DP(I) ,PITEV(I) 
GO TO 100 5 
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WRITE (4,624) JORKI (1,1) , FORKI (1, 2) Cel ШЕН 
НЕВІТЕ (4,625) ШОВКІ(І,1) ,ЛОБСІ (1,2) 29 (Т) ЕЛТШЕР 
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FLOW CHARTS 


І. MAIN. PROGRAM NO INPUT OR OUTPUT VARIABLES 


FILE DEFINITIONS 
ОЗ 


OG 
og 





READING THE TERMINAL 
WRITING To THE TERMINAL 


DISK STORAGE, DULT DATA 
DISK STORAGE, PERFORMANCE 


CALL 
NO 
1 Вар 


CALL 
COMPUTE 
"| COMPUT 3 





E CALL 
E > 










COMPUTE 


EXIT | 


STOP 
PROGRAM 





БАШ 





STOP 
PROGRAM 
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Z. BUILD SUBROUTINE 
THE RE ARE No INPUT OR OUTPUT VARIABLES 
FoR THis SUBROUTINE, HOWEVER SUBROUTINES 
CALLED BY THE BUILD SUBROUTINE DO 
HAN OLE INPUT AND OUTPUT DATA. 


CALL INST 
DETERMINES INSTRUCNONS 
DESIRED Br THE USER АМО 
TYPE OF TERMINAL USER | 


/5 OPERATING 


CALL. SYSTEMN 
THE SYSTEM CLASS/FIC ATION 
/S DETERMINED BY THIS 
SUBROUTINE. THE SYSTEM 
CLASSIFi CATION SHOWN 
IN FIGURE 2.6 15 


REPRESENTEO BY THE 


INTEGER VARIABLE, CLASS 
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SET АР IDENTIFICATION NUMBERS 


SELECTED 






FOR SYSTEM 
ссср 274 6 ©), CLAS 


CLASS T 12 
CESS E = 
GRASSES E 
CLASS ane g 
CASS V 2 
CUASS YL 5, 


SEE THE PRELIMINARY SECTION OF THE 
USERS MANUAL СОЮ EXPLINATION OF 


IDENTIFICATION NUMBERS. 






СЕОМ (1) = 113101 
G€oM(2) * 124 00| 
GEOM(3) = 136101 
GEom(4)= 14500} 
BRANCH ? 4 


&EOM (1) = 312201 
GE OM (2)= 323101 
GE om (3) = 324001 


GEOM (4): 335101 
GEOM (5)« 345001 
GEOM (6) : 356201 
BRANCH = 6 


авом (1) =512201 
вом (2): 5230! 
Geom (3)=525001 


Geom (4)5535101 
GEO (5)* 5562.01 
BRANCH = 5 








G&om (1): 2122.01 
GEom(2)= 223101 


GEom(3) - 224001l 


Geor (4): 236 0I 
СЕОМ (5) = 245001 
BRANCH = S 





GEOM (1) = 413 104 
Geom (2): 42400! 
Geom (3): 435101 
Ge0m(4)= 445001 
Geom (5)= 456201 
BRANCH = 5 


Geom (1) = 61310] 
GEOm (2): 625101 
GEOm(3)= 635101 


Geom (4): 6562.01 
BRANCH = 4 








CONTINUE 






DO FROM 





I=-1 TO 
BRANCH 


CALL MENU (M, TERM, TYPE, GCOM (I) 
INPUT: M, TERM, GEOM(r) 

OUTPUT: TYPE 

THIS SUBROUTINE DISPLAYS A MENU TO 
THE USER. THE USER MAKES A 
SELECTION FROM THIRTY FITTING 
CHOICES. THE FITTING ID NUMBER 
FOR THE NEXT FITTING TO BE INPUT. 


YES 












YES 
WRITE WARNING 
TO USER INPUT 
/ NO 
NEW ID NUMBER 
М=М+ | 
|5 
МО 1. BRANCH 
f 
YES 
K 






CALL SELECT (M, SORL, GEOM (I), TYPE, WORKI, WORKR) 
INPUT M, SORL, GEOM(L), TYPE 

OUTPUT: WORKI , WORKR 

THIS SUBROUTINE BRANCHES To THE FITTING 

SELECTED FROM THE MENU. 






INFORMATION 






ABOUT THE FITTING 15 CONTAINED ІМ 


THE ARRAYS WORKI Амо WORKR 


дом 1) = GEOM (L) + | 


CALL. SuMOu T (WORKI, WORKR ,M) 
INPUT: WORKI, WORKR, M 
OUTPUT. NONE 
THIS SUBROUTINE WRITES A Duc РАТА ЕБЕ 
ON THE USERS DISK. 


RETURN TO THE MAIN PROGRAM 
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JI. EDITING SUBROUTINE (ED) 
ITE IKE TSO INPUT OR ОатРит VASIABLES 
O CE e /( ОИТ МЕ HOWEVER SubRouTwes 
CALLED BY THE ED SUBROUTINE OO 
HANDLE INPUT AND OUTPUT DATA. 


READ THE 


DUCT DATA 
FILE 








ASk THE USER IF IT IS 
DESIRED TO CHANGE, DELETE, 
OR ADD TO THE DUCT DATA FILE. 


ПЕЩаНАНЕЕ e 


(ee Pele re С 


OUTFUT EDITED FILE 
CELL сам бит 


RETURN TO MAIN PROGRAM. 





WHAT 15 THE INDEX NUMBER, M, 
OF THE FITTING IN THE DUCT 


DATA FILE TO BE CHANGED 


USING THE MENU IF NECESSARY, 
SELECT THE NEW FITTING DESIRED 
AND ENTER INFORMATION 50 THAT 
VARIABLES IN Duct DATA ARRAYS 
With INDEX SELECTED PTS 
BE CHANGED. 


CALL MENU 
CALL SELECT 





MORE VES (в) 
CHANGES 
7 
К0 
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WHAT IS THE INDEX МИТЕЕК M. 
OF THE FITTING TO BE DELETED? 












FITTING TO SE 
DELETED THE LAST 
FITTING OF THE 










REWRITE THe 


FILE OMATTING 


THE LAST 
FITTING 








KEWoRK THE FITTING LD 
NUMBERS WITHIN THE 
AFFECTED BRANCH SO 
THAT ТО NUMBERS ARE 
CONSECUTIVE WITHIN THE 
BRANCH AND ARE ASSOCIATED 
WITH THE PROPER DATA AND 
\$ 










UNWANTED FITTING 
DELETED. 









WHAT IS THE INDEX NUMBER, M, 
OF THE FITTING THAT THE 
FITTING IS TO BE ADDED 
AFTER’ 


OPEN UP THE DATA ARRAYS TO ACCEPT THE 
NEW FITTINGS DATATA ГІР СЕРІК 

LOCATION, INDEX Mr Е а 
AFFECTED TO NUMBERS IN THE BRANCH 
WHERE THE FITTING IS ADDED 


USING THE MENU IF NECESSARY, INSERT 
THE DATA FoR THE NEW FITTING AT 
THE M+! INDEX 

CALL MENU 

CALL ЗЕ ВЕСИ 









ADD ANOTHER 


© 


TW. COMPUTE SUBROUTINE 


READ THE 
DUCT DATA 
FILE 


FIND OUT TRE SYSTEM CLASSIFICATION 
NUMBER 
CLASS = WORK I (1,1) /100 000 





BASED ON THE SYSTEM CLASSIFICA'T (ON 
NUMBER, SET UP THE APPROPRIATE 
NUMBER OF BRANCHES 









TO MATCH THE SYSTEM EATRACT 


FRom THE Duct DATA THE AREAS 
OF CONVERGING AND DIVERGING 
WYES (ADWB, ADWC, ADWM, ACWB, 
ACWC, ACWM) AND ANGLES OF 

FUNCTIONS (ALFAC,ALFAD) 







LOCATE THE ARRAY INDEX NUMBERS 


OF THE DUCT DATA FILE CORRESPONDING 
To THE FIRST FITTING OF A BRANCH 
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GET THe AMBIENT OPERATING CONDITION AND 
DESIRED POWER SETTING OF THE ENGINE 


(В) CALL PWRPT 


A FAN IS NOT INSTALLED 
IN CLASS 5 OR G SYSTEM ii 


<> YES 


GET FAN 
CHARACTERISTICS 
CALL FAN 











GO TO THE SYSTEM ШЕЕ Ne 
TAILORED TÒ COMBUTE ЕСЕ SE 
LOSSES FoR THE САИ МЕ 
BEING EXAMINED. 






ра 
сл 
|- 


р ОЕ 


SYST 
аа» CAE 
55/2 
CLASS 3 CALL 
SYS 3 







YOU WANT 


TO COMPUTE WITH 


A DIFFERENT 


POWER 
POINT 
үч 


NO 





KETURN TO 
TME MAIN 
PROGRAM 





فم 
сл‏ 
го‏ 





CALL 
BS 


NIS 
5155 


CALL 
SES 










YES 


ЕМА 5574 


Y SYSTEM THREE MATCHING SUEUR sS 
THIS SUBROUTINE \S CALLES BY TUE COS 
SECTION OF THE PROGRAM. ALL VARIABLES 
ARE INPUT FROM COMP SUBROUTINE, THERE 
ARE NO OuTPuT VARIABLES RETURNED TO 
СОМЕ? ALL OUTPUT 15 WRITTEN ORTIS 
PERFORMANCE FILE. 


SET UP Tue STARTING ANO STOPPING 


INDEXES FOR THE DATA ARRAYS FOR 


THE BRANCHES 





INITIALIZE SYSTEM VARIABLES FOR 
START OF 






ITERATION 











DUCT LOSSES 
INLOSS = 4.0 CINCH WG) 
Extoss = 6.0 (INCH we) 






Ерастое GAIN 
GAIN = -30.0 (PSF) 














COOLING FLOW РАБЪАСЕ 55-5 
Loss 500-052) 
COOLING FLOW 
Wee CEMMAX * RHOSTD /60.0 
BRANCH INFORMA TION 
рр45 = 100.0 (PSF) 
DP 56 =100.0 (PSF) 
тмоь =710 о СЕ) 
PTS = Рф ж 144.0 + DP56 





F3 







WASTE HEAT 


YES 
BOILER INSTALLE 


? 


GET ENGINE PERFORMANCE PARAMETERS 
BASED ON AMBIENT CONDITIONS, 


DESIRED HORSEPOWER AND RPM, AND 


INITIALIZED OR CURRENT FLOW 
ITERATION 


CON DITIONS 
SYSTEM. 


OF THE 





CALL ENGINE 


рэ 
Сл 
Е 


INITIALIZE INLET CONDITIONS FOR 
p RANCH N 

Воо 

РМА іу = РТ5 +1055 

PSEC. = ЕБУ САИ 

ртім = Рф х 144, 0 

ТШ КО И 


BRANCH 1-72 HAS COMBINED COOLING 
ANO ENGINE A\RELOW 
М = МС + МС 


FOR EACH FITTING IN BRANCH 1-2, 
CALE Е | г РС 
ТО COMPUTE PRESSURE LOSS IM 
THE FITTING. THE SUM 
DPIZ2= DPIz2 *DELP 
IS THE PRESSURE LOSS FOR TUE 
BRANCH 













NODE 2 IS A JUNCTION. COMPUTE 
TOTAL PRESSURE AND DENSITY AT NODE. 
PT2 = Pé x144 O — DP12 

RH02 = (PT2 - PV)/ K Tó+454.7) x R.) 





COMPUTE THE AVERAGE VELOCITIES 






IN THE THREE 
BRANCHES ENTERING AND LEAVING NOOE 2, A DIVERGENT 










WYE. 
BRANCA COOLING AIR: VDWB = WC /C Ru02 x A DWB) 
COMBINEO INLET: — VDWC s(wc«w2)/ RuoZ x ADwC) 

VOWm=: wz2/(RHOZ * ADWM 











MAIN ENGINE : 





COMPUTE NODE 5 PARAMETERS. NOOE S (S A 
CONVERGENT WYE, MIXING STREAMS OF DIFFERENT 
TEMPERATURES. IF NO WASTE HEAT Зо(в 15 
INSTALLED TEMPERATURE OF THE MAIN BRANCH , 






EXWAuST FROM THE ENGINE 15" 
EMANE TE ELSE, 
TMAIN = 170.0+(370 * 1072 x HP) 





COMPUTE TEMPERATURE IN COMBINED EXHAUST STACK 
BASED ON MIXING ENTHALOPY OF COOLING AND 
EXHAUST STREAMS. 
COOLING ENTHALOPY: HSEC=(1.421385E-Sx TMOD + 
221041) & TM0D + 5.6373 
EXHAUST ENTHALOPY: HMAIN * (.560SIE- 5 * TMAIN + 
22388) ХТМАІМ + 475396 
COMBINED ENTHALOPY: HSTACK = (W8/(W8-WC)) * 
HMAIN - ( WC/ONB *wC)) ж НЕС 
EXHAUST TEMPERATURE: T5 = (0.000841) # H STACK + 


4.33577} х HSTACK- 9.5778 


1S6 







COMPUTE ASSUMED DENSITIES AT NODE 5 
R40CC : PTS /(R*TS) 
RHOCB = PSEC /( R* Tmop) 

RHOCM * PmAIN / CR TMAIN) 






COMPUTE VELOLITIES AT NODE 5 BASED 





ON ASSUMED OR CURRENT ITERATION DENSITY 
VCWB = WC/CRHOCB * ACwB) . 

Це сет на +8) / (Кносс X ACWC) 

VCWMs W8/(RuOCM * ACWM) 






CoMPuTE THE VELOCITY PRESSURES AROUND 
NODE 5. 
PVB =RHOCB XVCWB x 2/(2.0 * 32.114) 
PVC: RHOCC* VCWCxx2./ (2 0 & 32. 114) 
PVM-RHOCMX VCWM «x2 /(2.0 x 32.174) 


COMPUTE THE STATIC PRESSURES AROUND 
NODE 5. 


РВ TSEC TD 
PSC: ENS EVE 


PSM = РМАМ- РУМ 





COMPUTE DENSITIES AROUND NOOE 5 


bASEO ON STATIC PRESSURES 
RHOCBT= ?58/(Вж Тмос) 
КОО ЦЕ 15) 
RHOCMT: PSM/CRX TMAIN) 









BETWEEN THE ASSuMED 





NO 
DENSITIES AND СОМРИТЕР Е 


VEN SITIES LESS THAS 


LET ASSUMED 
DEN SITIES EQUAL 


COMPUTED DENSITIES 
YES 
INITIALIZE INLET CONDITIONS E 


FOR BRANCH 5-6 
TIN: TS 
W=WC+W8 
РТМ = Рф ж 144.0 + DPS 
DPSG=00 


FOR EACH FITTING IN BRANCH 5-6, CALL FIT DP 
To COMPUTE PRESSURE LOSS IN THE FITTING. 





(Ез 
сл 
OO 


THE SUM рР56 = DPSE+ DELP 15 
THE PRESSURE LOSS FOR THE 
BRANCH. 











EXIT PRESSURE FROM BRANCH 5-6 
MUST BE ATMOSPHERIC, A BOUNDARY 
CONDITION COMPARE COMPUTED 

OUTLET PRESSURE WITH ATMOSPHERIC 
TEST = ABS(PTOUT- PD x 144 0) 








NO 


NEW INLET PRESSURE 
YES IS COMPUTED WITH 
LAST DCSE 











COMPUTE TOTAL PRESSURE 
AT NODE S 
PTS = PØX 144 0 + DP56 





INITIALIZE INLET CONDITIONS 
FOR BRANCH 2-3 

РТ АБ. 

TIN: TO *454 7 

ОР23-00 


FOR EACH FITTING ІМ BRANCH 2-3 
CALL FIT Dp TO COMPUTE PRESSURE 
LOSS JN THE FITTING. THE SUM 


DP23 = DP23+DELP 
IS THE PRESSURE LOSS FOR THE 
BRANCH 


INITIALIZE INLET CONDITIONS For 


BRANCH 3-5 
PTIN = P8 *¥ 144.9 


TN- T2 
DP35- 0.0 








FOR EACH FITTING IN BRANCH 3-5 
CALL FITDP TO COMPUTE PRESSURE 
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THE SuM DP35 = DP3S + DELP 
IS THE PRESSURE LOSS FOR THE 
BRANCH 









COMPARE EXIT PRESSURE FROM BRANCH 35 
WITH TOTAL PRESSURE AT NODE 5, PT 5, 
COMPUTED EARLIER 

TEST - ABS (PTOUT- PTS) 


NO 


NEW INLET PRESSURE 
IS COMPUTED WITH 
LASTS 


SINCE ASSUMED INLET AND EXHAUST 


LOSSES WERE USED TO COM PUTE ENGINE 
PERFORMANCE, COMPARE ASSUMED AND 


COMPUTED LOSSES AND REPEAT ITERATION 
IF REQUIRED 





TEST1 =ABS (ppi2 + DP23 - INLOSS) 
TEST2 =ABS( P35 + DPS 6-EXLOSS) 







DESEE RESTER YES 


THAN 1 PSF 








USING Com PUTED 
INLOSS AND EXLOSS 
REPEAT ITERATION 





NO 


INITIALIZE INLET 
CONDITIONS FOR 
BRANCH 2-4 


ВТ = WEZ 
ОР 242 с 00 
TIN = ТФ 





POK EACH FITTING IN BRANCH. Z-89, CALL 
Fir oP TO COMPUTE PRESSURE LOSS IN THE 
FITTING. THE SUM DP24 = DP24 + DELP 

IS THE PRESSURE LOSS FOR THE BRANCH 
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INITIALIZE INCET CONDITIONS 
FOR BRANCH 4-5. INLET 
PRESSURE 1S A FUNCTION OF 
FAN CHARACTERISTICS AND 
FLOW 
INS 
PTIN = PTOUT +(K ж (( 
WC /RHOSTD * 60.0)- 
СЕМ МАХ) ХХ 2 + 
DPMAX) 5.14646 
DP45 = 0.0 


FOR EACH FITTING IN BRANCH 
4-5 CALL FIT DP TO COMPUTE 
PRESSURE LOSS IN THE FITTING. 
THE SUM ОР45 - DP45 +DELP 
iS THE PRESSURE LOSS FOR 
THE BRANCH 












EXIT PRESSURE FROM BRANCH 
4.5 MUST BE EQUAL TO РТ5 
TEST = ABS (PTOUT -PTS) 








NO 


YES | INCREASE THE ASSUMED 
COOLING FLOW A SMALL 


AMOUNT AND REPEAT 
ITERATION 
CWC +O 





SYSTEM | MATCHED 
CALL OUTPUT AND PRINT 
SYSTEM RESULTS То THE 
USERS DISK 


RETURN TO MAIN 





PROGRAM 
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FITTING PRESSURE LOSS САВВА 
SUBROUTINE. SET uP TO COMPUTE 
PRESSURE LOSS AND VELOCITY DATA FOR 
20 FITTINGS LISTED IN THE MENU 


PA A 
GAS CONSTANT, R= 53-3424 
ROUGHNESS Е 1 0 000/5 


COMPUTE VISCOSITY OF AIR, MU, 
AS A FUNCTION OF TEMPERATURE 

маи = (29.95 +0./89304 юж TIN — 
2.67587 E-5 ж тум » 2) * 1-7 









ASSUME A VALUE FOR RHO 


RHO> PTIN/(R¥ TIN) 


COMPUTE VELOC OE ЕШ 


V* W/CRHO X DATA?) 


165 


В поло ЛИЙ д, ума | - 


COMPUTE PRESSURE VELOCITY 
PV = RHO xVxx2/(2-0 х 32.174) 
COMPUTE STATIC PRESSURE 
PS = PTIN - PV 


COMPUTE DENSITY OF GAS ENTERING 


FITTING USING PERFECT GAS LAW 
енот = PS/CRX TIN) 






COMPARE COMPUTED RHO AND 
ASSUMED RHO 


TEST = ABS(RHOT - RHO) 







TEST GREATER YES 


THAN 0.00] 


USING LATEST VALUE 
FOR HO COMPUTE 
VELOCITY OE FLOW 


COMPUTE REYNOLDS 


NUM BER 
ЕМ = ЕНОТ X DATA2 X V/Mu 
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GO TO THE COMPUTATION FoR THE TYPE OF FITTING 


IMPLE FITTING 
Е СОИ, 2827 


STRAIGHT DUCT 
# 2 


ELBOW 
#F 5,0,10 11,12 


SMOOTH RAD! RECTANGULAR ELBOW 
#7 


BRANCH SECTION DIV. WYE 
+ 13 


MAIN SECTION DIVERGENT WYE 
+ 14 


BRANCH SECTION CONV. WYE 
+15 


MAIN SECTION CONVERGENT WYE 
+16 


DIFFUSER 
#17, 18,19, 20 


CONTRACTIO 
ао 


INLET FILTER 
+ 24 


GAS TURBINE MODULE 
#26 


WASTE HEAT BOILER 
+21 
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SIMPLE FITTING, COEFFICIENT (DATA2) Times 






VELOCITY PRESSURE 
DP- DATA2 x pv 
PTOUT = PTIN - DP 
TOUT =TIN 








STRAIGHT DUCT, FRICTION FACTOR F, TIMES = 
TIMES VELOCITY PRESSURE L: DATAZ D:DATAZ 
F=0.25/(ALOGIO(CE/3.7 * DATA2)) + 
| 5.14/RN Xx X O3) * X2 
| DP = F X*DATA3/DATA2 ж РУ 
TOUT = TIN 


PTOUT= PTIN- DP 






ELBOW WITH REYNOLOS NUMBER 
CORRECTION FACTOR, KRE 


KRE -1.3101- 0.1485 *ALOG (RN x IE-4) 
DP = DATA3 x KEE x PV 


TOUT = TAN 
PTOUT=PTIN-pe 








i58 












SMOOTH RADIUS RECTANGULAR ELBOW 
WITHOUT VANES, REYNOLDS NUMBER CORRECTION 
FROM TABLE iN THIS SUBROUTINE 
RN< 10000 —* KRE = (Л 
Км > 10000 —> KRE FROM TABLE 
С = РАТА > 
рр: СХ ККЕ ж ру 
TOUT =TIN 
PY OuT = PTIN = DP 





BRANCH SECTION OF A DIVERGING WYE. LOSS 

IS DEPENDENT ON VELOCITIES IN MAIN SECTION 

(уруум), COMBINED SECTION (vpwc), BRANCH 

SECTION (VOWB), AND DIVERGENCE ANGLE 
TEST = (ADWM-Apwe 


K2 = (((-3.59147E-8 * ALFAD + аўсе 
3 8309E-6)X ALFA D+ 5.74 E-5) RR=VOWB/Vowm 


X ALFAD 4 |.04Е-3)х ACFA D IF RR>0.8 к1- 0.4 


4 0.000017 C= K1 «(10+ (VDWB/VDWC) xx 2- 
C Vo «(VbwB/vpwC)xx 2 - 2.0 2.0* (VbwB/vbwC)x Cos(AcrAp/ 
* (vow B/VDWC) * COS(ACFAI/ 57.3)) 
57.3)- KZ *(VOwW 9/ vpwC xx Z 





PV = RHO XVOWC ¥%2/(2.0 % 32-174) 
ОР = Ся РУ 

PV =RHO журмВжх2/(2..0 ж 32.174) 
PTOuT = PTIN- РР 

TOUT =TIN 




















MAIN SECTION OF A DIVERGING wre 
COEFFICIENT BASED ON VELOCITY RATIO 
C= 0.4 x« (1.0-VDWM /VDWC) xx Z 
PV = Кнож Урмсжж2 /(20жзо 114) 
DP = C РУ 
PV= RHO XVDWM ¥*72/(2.0% 32.174) 
ProuT « PTIN - DP 
TOUT UN 





BRANCH SECTION OF A CONVERGING WYE 
С.о t (VCwB/VCwC) Xx 2-20x ACWM 
ACWC X (VCWM/VCWC)XX2 - 2.0 x ACWB/ 
ACWC * (VCWS/VUCwC)*Xx2. x COS(AcFA C/ 53:1) 
РУ 2 RHOCC X VCWC X x 2/(2:0 X 32.074) 
DP s C * PV 
PTOUT = PTIN - DP 
TOuT = PTrouT/CR*RHOCC) 





MAIN SECTION OF A CONVERGING WYE 
KI= I. I$ X ( ACWB/ACWC)X»x* 2.21 + (1- 
ALFAC/60.0) X (ACWB/ACWC) X 0.4 
TEST = ABS(C ACWM- ACWC)/ACWC) 
(С TEST < ОО ПЕ К, - Оо 
С: |.о + (уСи М/усмес) * X2-2.0 X ACWM / 
ACwC X (vCwM/vCwcC)* 2.-2.0 ж 
(Acw В/А Сес) Х (УСмв /усис)Жх?2 ж 
С0$ (АсЕАС /57.3) Ec 
DP- C x PV 
Рт бит - PTIN- pP 
MOUSE TOU Ce occ) 
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DIFFUSERS 
1 =SQRT (4.0% DATA\ /3- 1416) 
КМ = RHO¥L*V/Mu 
LANDA з 0.3154/ ҚМ жж0.2-5 
СЕК = DATAZ X LAMDA 
С = СЕК + РАГА З 

pP = C x PV 

PTOuT = PTIN- DP 

TOUT =TIN 



















CONTRACTIONS 
V=W/CRHO * DATA4) 

Ру = Кножужж?/(20ж32.174) 

DP = РАТАЗ * PV 

PTOUT = PTIN— DP 

TOUT =TIN 








INLET FILTER, ECE Ее 
BASED ON FACE VELoc ITY 





PTouT = PTIN- DP 
TOUT = TIN 








рэ 
“| 
кз 





DP -(DATA2 X VX X DATA2) X 5.14696 



























GAS TURBINE MODULE , 6t LM?2500 
LOSS BASEO ON MASS FLOW RATE OF 


COOLING AIR THROUGH THE MODULE 
DP= (\.61E-3 XW **2.15) * 5.14646 
Py = 0.0 
PT OUT = PTIN - DP 

ОСЕ ИМЕ ТУЫ ү) x (He/20000.0 X 100.0) 


WASTE HEAT BOILER, GAS SIDE PRESSURE 
LOSS BASED ON TUBE BUNDLE GEOMETRY, 
VELOCITY JN NARROW PASSAGE OF TUBE 
BUNDLE, FRICTION FACTOR „F, ( FUNCTION 
OF REYNOLOS NUMBER 

V= Ух РАТА 7 

Pv= RHOXV жж2-/(2.о ж 32.174) 

KN = RHO x V x DATAS /mu 

F = 1.334 x EN % x (-0. 1453) 


DP = FX DATAZ X PV 
Ргоит = PT N - DP 
ТоЧТ = О. 00085 X HP 247.0 + ТФ 
RHO = (PS- DP) /CRx TouT) 

V* W(RHO x DATA |) 

Ру- RHOX VXXZ / (2.0 X 32.U19) 


RETURN TO SYSTEM! 
SUBROUTINE 
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USER'S MANUAL 


А. GENERAL 


The purpose of this program iS to analyze a marine yas 
turbine installation on board a ship complete vill ms 
exhaust, and cooling ductwork. Tae duct yeometry must ре 
input to the program to accomplish this. Ihe program makes 
a file called "duct data" which contains resistance informa- 
tion on each fitting entered. This file may be edited with 
tne built in editor cr if the user is satisfied with the 
Current design the file is read by tne program and used in 
the COMPUTE section otethe " ресарзаше COMPUTE uses tne duct 
data file and inputs dealing with the operating point of the 
engine to produce the performance parameters of the system. 
Performance includes Loth engine parameters and duct losses. 
All procedures in the program are accomplished using ап 
interactive terminal session. 

There are two versions o£ the program discussed in this 
user's manual. Version 1.0 is implemented on the NPS IBM 
3033 computer. Version 1.1 is implementel on the NPS VAX-11 
computer. 

This user's manual will discuss the questions posed by 
the rrogram. Faniliarity with the prograr sections and the 
questions asked in each section will facilitate program 
execution and help produce reasonable results. The most 
critical area for familiarity is iû the BUILD апа m 
sections of the program. Tt is not so critical аа 
COMPUTE section of the program because only two questions 
are asked for each operating point run after the aukient 


conditions are input. 
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B. PRELIMINARY 


The program does not design ducts or read mechanical 
drawings. The user lays a vital role by interpreting the 
system zcr the program. Some fittings are easy to recognize 
such as é€lbows, straight duct,  transistions, diffusers and 
contractions. Some are harder to understand, like diverging 
and converging  wyes. Each ritting listed in the menu is 
sketched for the user. The sketches show a typical view tut 
remember that the dimensions shcwn on the drawings are vari- 
able inputs so the configuration can change drastically by 
looking at a fitting over the гапче of variable dimensions. 

Before running the program the user should become 
familiar with the fitting sketches. Comparing the sketch to 
the fitting to be modeled will assist the user in preparing 
a list of fittings fcr the system. The user should note the 
dimensions and be prepared to input them to tne progran. 

The program looks for fittings in a definite sequence. 
Branches are groups cf fittings or sections of the ductwork. 
Branches run from node to node. A node is an entry, exit, 
ponction, (ап, Or engine. Refer to figure 2.6 for the 
various system configuratiors. Nodes are indicated in this 
figure Ey tne numbered black dots. Nodes have numbers from 
one to six. The branches get their number designation Егоп 
the end point nodes. The user shouli become familiar with 
the system schematics then it will be easy to understand the 
order that the program will Бе asking for fittings: 
Branches are entered in a sequence from the lowest number 
node to the next lowest number node untii all fittings are 
entered. For example, a class three system enters branches 
men oll owing order; 1-2, 2-3, 2-4, 3-5, 4-5, 5-6. To 
assist the user when entering fittings the program displays 
the current fitting identification number on the screen with 


the menu. The ID number is a six digit number where the 
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first digit is the system class, the next two digits ac Е 
Franch number and the last two numbers are the sequence 
humber of the fitting in the branch. A terminal session has 
teen recorded and the printout annotated to skow tnis 
number. 

It would be helpful to pencil in the node nuabers in the 


system drawings. The following table may help. 


TABLE II 


Node Designations 


1 Main air inlet (engine only or combined) 
2 Cooling air inlet or divergent wye 

off maln inlet 
3 The engine 

À fan 


Cooling air exit or convergent wye 
with main exhaust 


6 Main exhaust (engine only or combined) 


The user should prepare a list of fittings organized by 
tranches and continuous with regard to the sequence of 
fittings- It's the old "toe bone connected to the fcot 
tone" idea. As an example, the following list nay help. 

node 1 
vert intake, 3 orifices, with louvers 
straight duct 
rectangular contraction 
smooth radius rect elbow 
node 3 
etc. 
Do not forget to include abrupt exits where they aprear. 
Sometimes it is easy to overlook an obvıous fitting such as 


the engine module as part of the cooling air ductwork. 
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Only the class one system does not have either a diver- 
gent чуе or a convergent wye. class three and five have 
both. The divergent wye 1S fairly straight forward. The 
user only needs to enter the areas indicated in the sketch 
and the angle of divergence (0-90). The branch section of 
the divergent цуе із the first fitting in branch 2-4 (2-5 if 
no fan) andthe main section (combustion air) is the first 
fitting in branch 2-3. Ihe combined area and the diveryence 
angle are data entered when entering the branch of the 
diverging wye. The convergent wye is a more complex. It is 
located at node five. The Eranch of a convergent wye should 
БЕЛЕ last fitting cf branch 4-5 (2-5 ir no fan). It will 
usually be the fitting after the module. Tne main section 
(engine exhaust) of the convergent wye is the last fitting 
su branch. 3-5. Usually there are just two fittings in 
ppanch 3-5. The first is the nozzle or extension bolted to 
the exhaust plane flange of the engine, and the last is the 
Main section əf the convergent муе. The combined area and 
convergence angle are data entered with the Егапсп section. 
The convergence angle is usually zero ana the combined area 


is abcut equal to the sum of the aain and branch areas. 
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FITTING NAME: | NUMBER 
Vertical intake shaft with side orifices, | 01 
with or witnout louvers 


в. 














Se 


INLET 


aj 


| 


d 


AREA OPENINGS 


| | | Bl Г 6 q 

WITHOUT LOUVERS 2 2 3 4 
a = 

n 





| 

| 

| 

SHAFT 
| 1 

| سے‎ 
| WITH LOUVERS 


| ее | 
D————D P————————— СЫЗ ыб et Dtnáüós e liLil1iL(Óo sai 








| INPUT REQUIREMENTS: 7 
1. Тһе nuuber of orifices (152,5 9 TN) 

| 2. The cross section area of the Vertical Shae | 

3. with two orifices, whether they are adjacent 

cr cpposite | 

ц. If there are louvers installed 
| DUCT PATA FILE ENTRIES: — 
| ИСЕКЕ (Т, 1) WORKR (I,2) 4ORKR (I, 3) WORKR (I,H4) 
| Shaft 0.0 resistance shaft 

area Coefricient area 


rm ——— 


REMARKS; 
The lcuvers are flat plates of standari configuration. 
The opening areas are not required but should be 
ap rar mites; proportional to those Shown ІП БЕ 
sketch. 











| REFERENCE: = | | | 
Handbock or Hydraulic Resistance, Idel'chik 
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کک ست ےس 





Va — rra a sam — a t m RD ED TT TT s ee s s r AD ست‎ лы. 
| FITTING NAME 1 чомвкк: | 
Straight d 2 | 


SKEICH: 


D CIRCULAR 


RECTANGULAR 


| 


LOADS 
с 
i A / 
| қ 


1. Round: diameter and length . | 
2. rectangular: Е dimensions (a,b) 
eng 


| 
Б REQUIREMENTS: даки ра а. m | 





| DUCT DATA FILE ENTRIES: Б m 








| WORKR (I, 1) WORKR (I,2) WORKR (I,3) WORKR (I,4) | 

| section diameter or iength secticn | 

агеа equivalent area _| 
ВЕМАВК5: E NH H 


peace tes bach Equation used for resistance. . 
Friction factor ру correlation by Swamee & Jain. 
Equivalent circular diameter conputed for rectangular 
Sections. Length should be measured to the center of 


: - | 
Short fittings and to the start or end or a long 
fattig- 
REFERENCE: ` | in in 
Mechanics of Fluids, Shames 
E а. 
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FITTING NAME: | ОТ NUMBER: | 
Smooth radius round cross section elbow 03 | 


Ee SS ЈС 2 реси анне раны ламаны саамы ES ee Ee «чыр: 


ОКЕТСН 


Ss ee ць 


`N 


THETA 





LS аа аа  — P. „шш —— — a F — —  ——————— D 
H 
2 | 








PUT REQUIREMENTS: ` 

1. Cross section diameter f 

2. Radius of the turn measured to the centerline of 

the section | 
3% Ihe turn angle 
= се. —— ———— —— ——— ——— —1 
DUCT DATA FILE BN Wades: 

| ИСЕКЕ (Т, 1) 4 ORKR (1,2) НОВКЕ (Г, 3) ХОВКЕ (1,4) 
REMARKS: — ` J = _ M 


section ЕЕ Ее section 

area Coeficient area 

Turn angle should be from 0 to 90 degrees. | 
| 


a . ——— ——— 


= 
| " 
шыні І 
Nm | 
Ба 
ОН: 
uU 
LR 
гл 
“М 
са 
u^ 
сә 
= | 
‚< 
г] 
к 
rj 
ты 
r3 
(л 
о | 
ко 
со 
Z5 
+ 
Q 
ЕТ 
Ф 
м3 
ct 
(D 
"n 
w | 
“| 
| 
| 
| 


са, чш ш ы сш ES mas a A wm идите “mam O Е ЕСЕ 
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FITTING МАМЕ: NU 





| MBER: 
Segmented round cross section elbow 94 
3, 4, or 5 segments, 90 degree tura | 
KS E 
EKIETCH:; 


THREE SEGMENTS SHOWN 
(THERE MAY ALSO BE 
FOUR OR FIVE SEGMENTS) 


ee ey GS cee, Ge o — — o —— ФЕ — J تت‎ 


90 1 | 


D 


Numker of segments 

Cross section diameter | 
Radius of the turn measured to the centerline 
of the turn 


DUCT CATA FILE ENTRIES: P h N 
WOR KR (I, 1) WORKR (I,2) ПОХ КЕ (1,3) ГОККЕ (ГТ, 4) 
section 0.0 resistance Section 


area coefricient area 
EE — — ————- 
Note that the number of segments inciudes tne entry 


and exit segments. 


U! KEQUIREMENIS: 


— ты i ea سس ور سگ‎ i i а “am Sa аста" سے س‎ 
tt | 
= 
| 
x 


Е. аи 


б» 
г) 
Frjee 
c 
= 
с 
= 
ы. 
tj 
= 
н 
2 
t 
(л 
=A 
uo 
со 
=b 
4 
Q 
E 
w 
"O 
ст 
(D 
г 
CJ 
С) 
эы» —————— 


ا س ا 
morc eR oC‏ 
ÀSHE‏ | 
Mu 2-2 diua nan |‏ 
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| FITTING NAME: | | NUMBER: | 
Mitered round crcss section elbow 05 | 


| SKETCH: 


LL—————————E 





4 OPTIONAL CASCADED 
VANES 


THETA 


ӘЗ” тер ан 


ИТ РЕС ШТЕР ЕКЕНТ О 

. Cross section diameter 

+ IUCN angie | | 

. Whether or not concentric guide vanes are 
installed 





| 

| 

| 

| 

| 

| 

EUT FECUIREMENTS: о 
| 

| DUCT DATA FILE ENTRIES: | | Г | 
WCRKE(I,1) WORKR (1, 2) WORKa (I, 3) WOKRKR(I,4) 

| section diameter resistance Sectio | 
| ared coefficiert area 3 





` =p == Se GE حح‎ ee ee сЕ 





REMARKS: | S.C 
A Reynolds number correction is applied to this Zitting. 
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FITTING NAME: ! | NUMBER: 
Mitered rectangular cross section elbow 06 
without turning vanes | 


риша SEES OI DAMEN =“ <- сели Шо «ЕШ» ааш OO DU SS S ONU Sen ccc c  ) s ue СЕНЕ Еф cam ego ee ee ee дас сох дс л 


| SKETCH: 


| 
| 
W 
THETA 
| 





==  —  ———Á— — ras q r cnt 





eS trc mE a a с=з: с. EES “`ah твая? ше шы 





UT REQUIREMENTS: | 
- Height of the elbow, dimension parallel to turn 


axis 
. Width of the elbow, dimension in the turn plane 
. Turn angle 


ыы Пеле 








| DUCT DATA FILE ENTEIES: 
| WCRKE (I, 1) WORKR (1,2) WORKR (I, 3) НОЯКЕ (Т, 4) 
section hydraulic resistance Section 
area diareter coefficient area 
Macks C Snn ————— ше 


This fitting has a Reynolds nuaber correction aprlied 
to the resistance coefficient. 








DEFERENCE: 
| ASHRAE FUNDAMENTAIS 1981, chapter 33 


iom a AAA 26а нан ER AA все com ай 
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mm. m =з жый کس‎ чы „зз жы ОА тшшш ee — е. 


FITTING NAME; | [NN NUMBER: 
Smooth radius rectangular elbow witnout 07 
guide vanes 

SKETCHS 





= m ES. 








а a Чац P —————— w “sms. 








ааа 


bon БЕСОТВЕН ЕТ: ио 


1. Height of the elbow, the diuension parailel to the 
turn axis 


2. Width of the elbow, the dimension in the turn plane. 
3 Radius of the elbow measured to the centerline of 


| 
| 
| the elbow. 
| 4. Turn angle | 
DUCT DATA FILE ENT BEES: 
| WORKR(I,1) WORKR (I, 2) WORKR (І, 3) VORKR (I,U) 
section hydraulic resistance radius/ 
| агеа diareter coefficient width B 
REMARKS: ` E E | 
This fitting has a Reynolds number correction. 
The correction also varies with the R/W ratio. | 
а н — — 9 —— — ———————— А —— ——4 
REFERENCE: | 
ASHRAE ЕПНОАМЕНТАТ 5 1981, кБаргек 33 E 
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FITTING NAME: = | NUMBER: | 


Smooth radius rectangular elbow with 
Splitters 


(CEE canny MEE GS oe. ee EG ففق‎ RE 








ш» аана Sia Se Se —UPMEENEEO-CT i ИЕ qas > -—— n с 


ШЕШЕН: 


s= бл ie a ee ee ee ee ees Qe oe Ás 





TWO SPLITTERS SHOWN 
(THERE MAY ALSO BE 
ONE OR THREE) 











| INPUT EXQWIRERENTS:, |. 5) ar3 | 
. Number of splitters, ог 3 x 
| 





P 
1 | 1,2 
2 Height, dis ance parallel comturn axis 
3 Width, distance in turn plane | 
Ц Radius of elkow to section centerline 
5 Turn angle 
| DUCT DATA FILE ENTRIES: NU й қ 
i WCRKE(I,1) WORKR (I,2) WCRKR (I, 3) WORKR (1,4) 
section 0.0 resistance section 
area coefficient area | 
REMARKS: E ü К а | 
None 
1 
MEEPERENCE: о 2 d 
| ASHRAE FUNDAMENTALS 1981, chapter 33 | 
— c سس‎ | 
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FITTING NAME: | | NUMBER: 
Mitered rectangular elbcw with vanes 09 








eS س‎ U «ань flap ی و چ‎ 


SKETCH: 


— — —F. — ыыы 





THREE VANES SHOWN 


| 

(THERE MAY ALSO BE 
| ONE OR TWO) 
| 
| 
| 


| 
| 
1 
| 
| 
| 
x 


T" -————— AEE ee i салла» ч са BR Le, ® Ке = ee 


Number of vanes (1, 2, ог 3) 


UT BEQUIREMENTS:; 
я Cross section area 


| DUCT DATA FILE ENTRIES: 
| WORKE (I,1) WORKR I,2) WORKR (I, 3) КОБКЕ (1,4) 


Section resistance section 
area coefricient area 


қ)“. 
Е» 
р 
ct 
(0 
= 
с 
ч 
t3 
H. 
= 

те 
< 
fw 
r3 
Ф 
(n 
p 
ч 
m 
с 
(n 
(0 
fou 


SE — „ш з —— шшш шшш m= — шшш ар на шш й RR — EE EE ED 6 Quan sima “ара سسس سس سے‎ 


— s “ a u O G s “s a suntu D wal “ms <s ra t a m НЕ ucu s m Q aa rrr w cr 
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атым» нье ань Оет сьо аЬ ЦЬ льо оь НЬ 





{ FITTING NAME: = | | 104958: | 
Rectangular elbow with converging or 10 


diverging flew ii | 
| SKETCH: | 























| 
| | 
| 
| | 
| WO | 
| 
| | 
| | 
| | | | 
| М1 C. | 
| | 
INPUT REQUIREMENTS: 00 | p 
| 1. Inlet height, dimension parallel to turn axis 
2. pou Degmgnto dimension parallel to turn axis | 
| 3. Constant width, dimension in turn plane | 
DUCI LATA FILE ENTEFIES: ш га = 
WORXR (1,1) WORKE (I,2) WORKR (I, 3) КОВКЕ (1,4) 
inlet inlet hyd. resistance outlet 
area diameter coefficient area 
REMARKS: ee | 
Elbow should have a 90 deg turn. . 
The width should remain constant in the elbow. | 


ccÁ—— CE COSMAE EDI C ee ee ENNSD-UANEVE £30 c c ——"—-H"E— растанне» ананас стана EE ee Se ito oc mane Сі ста uds ec Se OS 


NCE: | 
AE FUNDAMENTALS 1981, chapter 33 | 


“£ 


-—— €————ÓÁ————— ла л — — l l ы ашады eee E MES reel 








I TING NAMES 


Two 90 degree rectangular elbows in a 


Z-shaped configuration 
ӘКЕШІІ: 


UT REQUIREMENTS: 
` Width of elbows 


DUCI ГАТ 
WOR KE ( 
Sectio 
area 


н 
` 





Height of elbows, dimension parallel to turn axis | 
dimension in turn axis 
The distance tetween the centerlines of the offset 


WORKR (I, 3) НОККЕ (Т,4) | 
resistance section l 
coerricient area | 


REMARKS: СҢ о 
Мопе. 


шт 
ша 
zm 
2» 
[т] 
кее 
c 
< 
о 
o 
< 
С] 
= 
lar 
о» 
L4 
п 





| 

| 

ee 
NCE 


198 1, chapter =н | 





—= анна: aa сабачага аа? س سے‎ ЕЕС 0 A 
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FITTING NAME: š . 
Two 90 degree elkows in different planes 





танг ИИН? ae ee RS ee CET Coma 20 ш сщу. 





SKETCH: 


Ac IE ERE Ы... 
С 
с“ 





ОЕ ОСЕТЕ ЕНИ: | | | 
Е а of elbow, dimension parallei to turn axis 
. Width of elbcw, dimension in the plane of the turn | 
. Distance between the centerlines of the duct 
connected to this arrangement 


Se SE ee nD ne ae ee ea ee ea ee ee ee Se точен со 


DUCT CATA FILE ENTRIES: 





| WORKR (I, 1) WORKR (I, 2) WORKE (I, 3) WORKRE (I, 4} 
section Hydraulic resistance section 
area Diameter coefficient area 

| ЕЕЕ e ол 


Resistance coefficient is a curve flit to the 
tabulated data. 


| 
L 





REFERENCE: 
ASHRAE FUNDAMENTALS 1981, chapter 33 | 


- 


қа заяви со шап ТО стои пра i—aaaaaiaaaku -ттеевацечато сее ст p ателие с ее a a 
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^ cem — = 











FITTING NAME: NUMBER: 
Diverging wye, branch amda пати сесте 13 & 14 
т ы... 

SKETCH: E 


اق رو ш‏ < 


COMBINED INLET AREA 
(ADWC) 





а 


BRANCH 
ы MAIN JET AREA ` 
MODEM (ADWM) 
INEUT FEQUIREMENTS: 3 М а 
А. Branch section В. Maın section 
1. combined area 1. mala area 


2.  Lranch area 
3. divergence angle 

















| DUCT DATA FILLE ENTRIES: быя Ue 
КСЕКЕ (Т. 1) WORKR (I,2) KORKR (IL, 3) WORKR (i,4) 
combine branch divergence рама С П 
агеа агеа angle area 
== RU TEN ананан ЕСЕ ----4 
| REMARKS; | 


Тһе divergence angle should follow some centerlire | 
streamline. The areas are cross section areas 
perpendicular to the streamiine in the sections 

away from the dividing location.  Cooiing,air flows 
through the oranca section. Main inlet air to the 
engine flows through the main section. Both ‘flow 
througa the combined section. 





— —  — — کد سے‎ A A-———M——M'J———— کے ھ ص‎ ааыа ЫЕ اک چ ست ن ا ف د‎ дыг 2 С. с әааиниБаныәашсдаадиР арәәлы: ӘӘ 


REFERENCE: | | : 
| Handbook of Hydraulic Resistance, Idel'chik 


E 


цц n e T. 20 
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FITTING NAME: | | NUHSER: | 
Convergent wye, Ггапсһ апа main sections | 15 6 1 















SKETCH: 
SCHEMATIC 
| MODULE EDUCTOR 
COMBINED 
AREA 







COMBINED 


BRANCH 
AREA 


(ACWB) 


tri 
ж 
ае 
> 
с. 
2 
3 


(RECTANGULAR 


MAIN CONTRACTION FOLLECTOR 


AREA 
(ACWM) 


a‏ ص 


атыша бй.‏ اا سے کاک ےک а, баі Фіна ge MIN a ES Со ОЩЕ ето Иса‏ ا ےو وک ہے مک ہے ی 





| 
— ы v — uror. у 
INPUI-REOQUIREMENTS: | | | 
А. Branch secticn Bey Marin section | 
12 on C n area 1. main area | 
| 2.  ccubined area 
3. | convergence angie 
| DUCI CATA FILE ENTRIES: fitting l5) 
| WORKE (I, 1) WORKR (1,2) WORKR (I,3) WORKR(I,4) | 
| сопріпе branch converjence branch | 
| агеа агеа angie area | 


| ranch area has module cooling air flowing through 
1 The main area has engine exhaust flowing ШЕШ 

| it. The combined area haS both. The angle Should fe 
measured to representative streamlines at the plane 

| where the two flows meet. 


| 


REF ERENCE: | | | 
Handbook of Hydraulic Resistance, Idel'chik 











— MÀ 


ee‏ > ي 
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FITTING NAME NAME: E = q NUMBEBE 
| Conical diffuser | 17 


l————M— — нання 
SKETCH: 


БӘРІ)! 


| 
| 
| 


та але аа Маа аа аа анна аа = p€————————————————— —— P s — чара اک‎ Ҹам  — im s Ge ae ae 


в 


(ләкә ән) 


ОТ БЕСОІВЕМЕМТ5: 

. Length of the diffuser 

ё Inlet diameter 

. Cutiet diameter 

š IS there distorted flow at the iniet 

. Аге tnere dividing wall ог bafrles installed to 
reduce resistarce 


2 





Sa SS SS Se a Ss em: Se me ee t ma as = a W .— a s . : fn 


| DUCT CATA FILE ENTSTES- 

| КСКАЕ (Т,1) КОЕКЕ (Т,2) WORE (1, 3) WORKR (CEP 
iniet friction LOW o outlet 
area coefficient Сое устеп area 








REMARKS: 
| The Rom recognizes a wide diverging angle and 
warns the user. Resistance in this сазе пау Бе 
reduced by 35 % with installation of baffles. 


| 

| 

= о а 
REFERENCE: | | 


| Handbook of Hydraulic Resistance, Idel'chik 


ےک ma... a...‏ سرد کت ہے س و 





ааа‏ = = == بے س ا 








— À HR = o —À M A аты еа 


FITTING NAME: Son 
Flane in-line diffuser 





SKEICH: 





SS SE SE EE VLL иван С: ama ee 


IG 
^u 


T REQUIREMENTS: | 
Length of the diffuser | 
The constant heigth of the diffuser 
The inlet width 
The outlet width 
Distorted flow 
Installation of baffles 





‘ANE CON) го 


e 8 е е е е Cj 


DUCT DATA FILE ENTRIES: 
WORKR (I, 1) WORKR (I, 2) WORKR (I, 3) 
inlet fr iction loa и 
area coefficient coefficient 


| REMARKS: 


The divergence is assumed to be uniform wita respect 
to the main centerlire. А wide divergence angle 15 
recognized and the user is asked if dividing walis or 
baffles are installed to reduce resistance. 


ШЕВЕБЕНСЕ: 


Handbook ог Hydraulic Resistance, Idel'chik 


a Se een cum. aum amps cts O Guns emm eee аа» сачы» «нан S, «ыны, ачы» «ва» «н» 
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ict 
- 


о i — É 





аьла ДА ьо ъа алар Quos anco Se 


| 


| FITTING NAME: е NUMBER: | 


Pyramidal in-line diffuser 





Serie ars 


за 








— 8 


UT FEQUIREMENTS: 2 | 
- Length of the diffuser | | | 

. Smaller inlet dimension, larger inlet dimension 

. Dimensions parallel to inlet dimensions 

. Non-uniform velocity prorile 

- Are baffles installed | 








= === rT 
j ЯСЕКЕ (1,1) WORKR (I,2) WORKR (I, 3) НЭБКЕ (1,4) | 
| inlet Irictios flow о. outlet | 

агеа coeificient coefficient area Е 
| REMARKS: SS ц ^ 


A unifcrm divergence with гезе to the centerline is 
assumed. Wide divergence angle is recognized by the 
POE With a wide angle the flow resitan c саш 

e reduced by 357 with bāffles or dividing wails. | 





| 
НЕН Тл LL 
REFERENCE. | 


| Handrook of Hydraulic Resistance, Idei'cnik | 


| E سے‎ ee el س سے‎ 





= 
FITTING NAME: | NUMBER: 
Transition diffuser, round to rectanguiar 20 
or rectangular to round 


SKETCH: 


BE || 


| | 
M 
x 
| p 
| | 
| | 
| 
| 
| 
| 
| 


REQUIREMENTS: 
Manner of tranSistion 
Diameter 
rectangular dimensions 
Length of the diffuser 
Non-uniform velocity distribution 
Installation of baffles or dividing walls 


T DATA FILE ENTRIES: 
WCEKR (1,1) WORER (1,2) WORKR (1, 3) WORKE (1, 4) 
inlet friction flow . outlet 
area coefficient coefficient area 


Э 
| 
| 
| 
| 


H 
= 


іе е е е е е с: 


EZ 











Е спо О 


REMARKS: 
Unifcra divergence with respect to the centerline is 
assumed. Wide divergence angle is recognized and if 
baffles oc dividing walls are installed the resistance 
is reduced by 35%. 





| REFERENCE: ` | | т 
| Hanirook of Hydraulic Resistance, Idel'chik 


— — —  — — — k  —— سے س‎ — Ul ÉL E. => x = = 


FITTING NAME: | 
Circular сопа е соп 








Length of the contraction 
Upstream diameter | 


UT REQUIREMENTS: 
. Downstream diameter ' | 








ТА FILE ENT HIES: | 
WORKR (I, 1) WORKR (I,2) WORKR (1,3) WORKR (1,4) 
outlet 0.0 resistance inlet 
area coefficient area 


If you,need a transitional contraction you could use 
this fitting or fitting 22. The area ol rp o pn > 
or outlet would have tO be converted to a Circle or 
rectangle as required by the geometry for input to 
the rrogram. 


l REMARKS: | 


кеі a ee CIE DC Le a A ee le ars s f naraU a <s 


REFERENCE. 
ASHRAE FUNDAHENTAIS 1981 ще паре 33 
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FITTING NAME: | 
Rectangular contraction 


Е Е Е. 
| ШІТБЕР: | 
| 22 


| 


men ИЕ ЕЕЕЕЕЕЕЕЕЕЕС ЕШЕЕЕЕССЕЕСЕЕСЕЕЕСЕЕССЕЕЕЕЕЕ С 
SKETCH: 


Bl L 


INPUT БЕСОТВЕМЕНТ5: | 
Т. Length of the contraction 
2. Upstream dimensions 
3. Downstream dimensions 








ный i, سے اق وک‎ . чш — ge 
| 
з 
p———————— J"— ———————————— ————— L AY T A S u uuu i 


— a asas ED ee ee a: ee — EPUM. ER Os FE ms a e 7 


DUCI PAPA FILE ENTRIES: 
| ИСЕКЕ (Т, 1) НОЕКЕ (Т,2) WORKE (i, 3) КОЕКЕ (Т, 4) | 
| outlet 0.0 resistance inlet | 
агеа сосы степь area | 


This fitting can be substituted for a transitioral. 
contraction. The inlet or outlet area should rezain 
the same and the area for transition converted as 

| required. | 





REMARKS; | 
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= ЖЕ =. т о ee ee WK sama s a n a sr UJ ll" 


FITTING NAME: 


Аа 
NUMBERS 
Screen 23 


SKETCH: 










SCREEN AREA 
(FREE FLOW MEANS HOLE SPACES) 





^ 


“DUCT AREA (OVERALL AREA) 


Аа. (абша Чач а ны. иран «о и MM «ава тра» Фіна. арыл, шьет. il ош 
pe————ÁÓ— ———————Á— Pb itd 


MÀ 


=] Se ee m Sm m = == —m. m ae a ee ee E m ma “= “sss Qam as мм 


INPUT REQUIREMENTS: 
. Overall duct cross section area 
2. Screen free flow area 





ТА МТЕІЕ5: 

ORKR (I, НОЕКК (І,2) WORKR (I, 3) WORKR (I, 4) | 

duct 2120, resistance duet 

area coefficient area 
REMARKS: 1 7 77 м ж 


This fitting is useful for the screen in front of 
the engine inlet. The free flow area is the sum of 
| all the holes in the screen. 





197 


OO Lh TF LTN 
FITTING NAME: NUMBER: | 

Louver inlet 24 

———Á  ————————————————————————— M ——— — 
ӘЖСЕТСН? 











DISTANCE 
ACROSS ` 
LOUVERS 


FLOW 


DISTANCE 
BETWEEN 
LOUVERS 





| 
| 
| 
| 
| 
| 
| 
INEUT FECUIREMENTS: ED m Š | 


Distance across tne louver openings 
Distance between the louvecs 


<— —  . — —— — a  — — =a . Á. L ——== = 
М) 
e е 


| 
⁄ | 
3. The number of openings between the louvers | 
4.  puct area just inside the louvers | 
DOCI LATA FILE ENTRIES: | 
НСЕКЕ (1,1) НОЕКЕ (1,2) WORKRB (I, 3) AORKR (I, U) | 

d UGE 0.0 resistance duct 
area coefficient area | 


REMARKS: | 
| Ihe correlation is for flat louvers with the front. | 
edges flat with the face of tne louvers. No friction | 
1s included in this correlation. setter models need 
to be Seed ge fcr louvers with moisture separator 
edges. The louver angle is 45 degrees to tne face. 
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Filt 


FITTING | NAME: o Î ugg: | 
| 25 








SKETCHS 





D a ЧЫЙ сь ———— husk‏ سے کک سے س ہے .و ےک ا سد کے ی 


| INPUT REQUIREMENTS; 

None if the default value is used. 

If another filter type is to be used then the 
user pO Е pressure loss data as a 
function of face velocity. Only a few points 
are required for the power curve fit to work. 
The number of points is an input (two min.) 

DUCI DATA MEME ENTBDESS 
WORKR (I НОЕКЕ (1,2) WORKR (І WORKR (1,4) 
ри ter P multiplier exponent filter face 





. 
ت 


vm 


| area (А ( B) area Я 


шь сщз сш ащ. —A———————————'L'a—————JR————— ee ee ee ae eee ee eee‏ س ا ا 


Ee z 
| REFERENCE: 
| Filter manufacturer's data 





pe a NN < ° “аа” eS GS odii: cd 
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FITTING NAME: a m " NUMBER: | 
Multi-Ltaffle type silencer 26 E 


ЕН: 


DISTANCE ACROSS DUCT 





x 
| 
| 
! 
| 
N TS НЕ NUMBER | 
OF GAPS | 
| | 
| 
bi 
INEUT REQUIREMENTS: | аа. 
Ша за between Lafiles 
2.  Baffle thickness, 
3.  Baf£rle length (with flow) 
4. Duct dimensicn parallel to gap | 
5. Cuct dimension across japs | 
| 6. The number of gaps а 
| eee. DATA FILE ENTRIES: | 
ИСЕКЕ (Т, 1) Ww ORKR (I, 2) WORKR (I, 3) WORKR (I, 4) 
duct 0.0 resistance duct 
area coefficient area | 














REMARKS; 
| This is a composit model. The resistance coefficient 
is medeled as a sudden contraction, friction along the 
length of the baffle, and a sudden expansion. it is 
not a very Еос mcdel and a model vased on experimental 
data woul e better. 








| REFERENCE: | | | | | 
| NAVSEA Inlet Design Handbook for Marine Gas turbines 


cere me ы БИЕ А ре ста‏ اا 


ССВ) 
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[ FITTING NAME: 
| Gas turbine module 


| ме 
| TM 


и MM MN 


SKETCH: 








GAS TURBINE MODULE 
x* COOLING AIR PASSAGES OSD IE S 


| 
| 
| 
| 
x 
| 


| 
| 
, 
44 








Е СЕИ REQUIREMENTS: 
None 


Г; 
| 
, 
| 
| 
| 
| 
| 
au 








| DUCT LATA FILE ENTRIES: 
ОЕТ WOBKE (1,2) ПО) WORKR (I, 4) 
| 1. 0 1.0 1.0 1.0 
| REMARKS: um й " i 


air through tne module. It is a power fit to data fren 
General Electric Co. It should be good as long as 
RO and exit areas remain about the same. The 1.0's 
in the duct data file are there to prevent 31 91 5100 Бу 
zero in the program and are not actually usei. 


| This model is based on the mass flow rate of cooling | 


REFERENCE; Пс 
Manufacturer's data 





> ————— P i ——— n — AAA vi Ж З , — < Ст سے کک سے‎ 
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SS == Е eee кес, 
FITTING NAME: . NUMBER: 
Waste heat recovery boiler 28 
SKETCH: 0000 Е ES: 





| 
| BUNDLE 
| EQUIVALENT 
| DIAMETER 
NUMBER OF мб. 
TUBES IN À —— KN 2 [DISTANCE BETWEEN 
ROW / A 7 ROWS 
t ——— — 
) \ / 


NUMBER OF ROWS TUBE | 
| 
| 

` P d ur | 


DISTANCE 
BETWEEN 
TUBES, INA 


| 

| 

| 

ROW 

| 

| | 
Наван ne T | 


A default is available. It is based on current design | 
| considerations in the racer program. Should you 
choose not to use it several inputs are required. 
Read the reference and be preparei to enter the values 
Shcwn on the sketch and a few you will have to compute 
| Wewmocl: (mee. Cube equiv. dia. and hydraulic dia.). 
DU 


ee ED. ee ee: ee ee | is Sm mm с шшш a ie ce 


ЕП DATA PILE ENTRIES: 











j WORKR(I,1) WORKR (I.2) WORKR (I, 3) WORKR (I,4) | 
| duct 0-0 resistance luct | 
агеа coefficient area | 
—— — — —— — —— —À ———— — — —— — —— ——— — —— o 4 
REMARKS: | i : | 
If the manufacturer will FE the data, write your 
own model, but this should be close for preliminary 
studies. | 
| 
= --------------------------------------------1 
| REFERENCE: | 
| Extended Surface Heat Transfer, Kerns and Krauss | 
| радес 562-589 | 
ААА а. аер ыыт у — 4 
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C. EXECUTING THE PRCGRAM 


1. IBM 3033 at NES 


Issue the following commands to compile and execute 

the program. 

FCRTEX filename 

СИВ АШ Ве воле те ро MOD2EEH NONIMSL 

ICAD filename (START 
"filename" is the name of the program in tne user's files. 
NONIMSL is required because the program calls the NONIMSL 
library with FRTCMS when defining files and clearing tne CRT 
screen. If the file has been compiled on the user's disk 
the lengthy compiling may be omitted and issue just the iast 


two lines. 


М а 


ler 
{2 


Р 


lin 


The program version to be used 1s 1.1. This version 
is a modified version or the program listed in Appendix A. 
The modifications include elimination of ail calls to 
ББС 5, FRICMS 15 used for two purposes in version 1.0. 
First to set up file definitions and second to clear the 
screen at appropriate times to prevent the format cf the 
display from being chopped up. Пре 31e deripitions іп 
version 1.1 are set up using the standard OFEN statement of 
FORTRAN 77 used on the VAX-11/780 at NPS. ALL аса рове 
FRTCMS to clear the screen were deleted ard are not needel 
on the VAX because it scrolls tne display гол the bottor 
and does not cut off any continuous screen displays. One 
other change was made in the file definition area, а 
writes tc the terminal where made to unit 5 and all reads 
from the terminal were made from unit 6. This agrees with 
the convention of FORTRAN 77 as implemented on the VAX. The 
рата сад runs like any otner program on the VAX, first the 


program must be compiled using the fortran compiler, then 
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linked and run. The program is stiil interactive on theme. 
and about the only word of caution required 15 to remember 
to use CAPS ON or upper case input гог Толіса: r ша 
Using lower case leaves the user in a ioop where the projram 
keeps asking гог for a Correct гер Ine duct geometry 
fale inrormatioa0 ison a file са! 1ед düct.lat ana ® 


perforaarce inzormation is on a fií.e са! ей ошеошестани 


D. BUILDING A DUCT CATA FILE 


The following pages are a recorded session at the 
terminal where the author entered a system in to tne 
SEO Idm The system modeled is made цр from drawings for 
tne propcsed Arleigh Burke class guided missile destroyer. 
ine session has been annotated to point out features of the 


program. 


205 





— 


заць ўс ON I 5: 


m 
іш 


¬ 


ЕЛДЕГІ 


4022 
PHEN d. 


FORT 40 D2 


сте 
.. 


BY LCDR. 






x 
` 
ur uj 
A 
2 со 
5% 
> 
< 
. tg е 
ка ° => x= (E3) 
а ж к Ut) е рда с 
x: Ж * g S Hire в, (9409 2) 
kl ра * ш 3 T лл е ө) 
Ен н с. с. w ¥ ZU Orb) m 04 
tn Uv) +3 = £4 q су =@ >” гэ {+ „л 
м ыы A ы ғ S š TE eH B) O 
n a н E 55) 4 м аж nie Баб. = 
ж () =, 2 umn тО н + O £4 
єз О d н 7-2 % E 27 HOZ bu - 
о f г. ч hase (50) filiu mu 
5 ry г. 21 > = | > >, ж oO га I) 
ea O Ou ~ 4 ы v “От тра (З tH 
> = “о. > «x осо > SIS танго Z [A 
Р] н н ` — H = > 2 < м Оба алт. еп. О 
чара 0 ы (n >; £d t) O m Zu NH О 
к U „3 ص‎ N m ËJ eum ЕО Он, 2410 H Еч 
H = “© - fx] ща = C O асоба Û 
(ou rx U uy “= «Ба E4 H . ж он O O гг: 
т min He о Uv) < я «б DEI ғаны UM >з 
зым MZ we =) 5264 ra > H Ор О м и T PE 
сыға O bit Fm Oli г. 1:1 21421 Бао ра Еч T7 
одат 4 ж Е О HH о fd м ЕЕ КЛ Мо о чо о 
pO dt > H H HH зе O WAE i тз THG б => 
ty ғи 235 HÀ E оа n fA 4 зона тт АА О — 14 
меча Qi) «4 ptr, Se wo се MO ж Ош Бао. п U 
абуха =a; > mul E = а фа но zm дон ж m m 
mem © OH t4 т; HA met H fx] ËJ он CIDE A044 I} 
LIC O ta O QU t4 бы нм O 4 Be она сурт) ЕР (ц > 
п OVU wa > м mE хт "7 - мы в сара м ж 
оо» Wo > him С” H a a Et помо MHU O t4 O 
ОЧЕН H H e hao оз E4 Е HHMU во HUY O 
«йг, (т; 0421 fx, ЗО кан H uU жи соев ем m О 
жн аш HE О t = ес ес moO HZ Enw ALIH «tH H aft 
ар: ZEH а ОҒ Е н H O WIX mer © е 
(uU f м O AM mu t5 e UO FE 9116431111103 ЕР! 
=E = - 4 oo OU OUZ = Luxe ар вис мар uro 
OWadvibiet 4) н © [du H H е кт заа Ор ера. Р Ore 
ed] Uam HH и. GE eth +4 Lt, Etu р не оно” BD = 
ены тэм m E о о f: HIJA г. EEO MEO 
ға = <f O = Z0 O O AOO 0:21 1Е4Е4ОЕ400 U мен 
а «в аа р) # £4 CE Qt) С) = ж шымғгыз оза [3 E4 
о Он OF < С H ен Наннан ЕЗ ес UU 
на DMO ¥ n = aO or E +3 - = ғооо G zn 
UNOCAL H > о әт I 2 о бо UfIDHoA;HUT MIE ANN 
Gi rH x pA FA О „1 WH HD = о шо ОО су пешжойао г.» 
HDA OZ A ` +4 cju uc) а о ғы неар ОО РІА ие 
tmfsu omo U Ert] E. bli шо О zZ eu WHA UNS TOU 
= т> м MNO on 5: =< pla О ТСЫЛТРИО ОУ” 
.. H P < к= «Н OZOOHtfüH H 
т .. = DT pe Іш is] Án оз Бч) И СаО а НЫМ mn (Л) 
>, а m O >59 hin p p; uyo; HOW kiret (OO тж; 
О O «а O D AU Imt] t! шч то] ра, Нн аа П HIG OP + O 
н т Œ но O жы д ж — ІП ж 3: ЗОО Ожичае tu UO CO UO 
Е ы пы s ле лын кі рї лс 
04 w ж па эц опы. fe 5 Ha 
O хі * HO о ом Oo n^ v OM NVON e аа ° е 
* LAOA O ц ме ма“ р“ м HH ыы" NMF м ص‎ 
a n1 о > > له‎ рвілімні 


ZO CONTINUES 


— 


GERO 


NTER 


D 
аф 
* 


> 











о а 
HZ = 
"о Гаў 
QU FH ЖО Ез 
рагач а U 
pira) аа = 
U WH zt > ша 
ITN 421,450) ~ 
оо ІІІ) ws m 
22 t >o, DH г 
жесі 2% sO e U 
+40: < нано т” = га с 
«пя: 1 «ЕЕ 6731 DI m bla m 
23 e MOM =E Gn) +) P) раз by 
мухи ие е Umi Hoó3 U* £4 
ма а 4 О ДО ОС) п н oH 4 
мк ІК ет «бес 1 02072 >, O „аи H 
hae EOE Mmm et D] ра # ра 
тж ға тт м Бра 
ын b) п. 0302004 6G HNO В 
(тт мо > КО НЕА 43 №40] 8-27 T 
гуру дра пара ғаға d (lD) Ф 4 
HOD CDI III MOU 09 => 5: 79 UP n) n 
со пил тъп оти) аэра Та 7551 fra tr: і г. 
17:0) 01:2 азота и FAO PI HIE: О H FA Uu 
Lt Оо, раа ЕЕ зра Б Еа ЕЕЕ (Z5 — à - 04 N 
Еа с а аг, vCTU)DQ."Je3 tdi 466065 «с сэ Uu) [£5 | fs) сез 
HOOH нччносоошоновамаа са f C4 «л п = C4 
(00003. QO]ODEUDUO-fsO0tEDpg G0 п са о t3 > Ез N 
uv “4 "x m (е; O = 
quot DAD rN TUONO FIN су O - .4 2 со 
eere reN NNNNNNNNNMMT га "a с. е Uu о г. 
UA U fr} - ^ m fd ~ A 
nnns o 398 4 948 4*8 A et н = ^ ofr) п т ta > ^ 
t) ` г. 4 ` fac) > FA Не N 
ш т; mm ° H > na Us = ва > m 
г за =€ fr 4 fu ~ ш 2..4 О = — 14 
а а A + mua Сей Р m AH J ni m m 
t4 =: = UV BH 52 PIH ол C5 >: aiu) n О o x 
9 uu 5 V ш>; © x ÉD £40 04 = D ын m ге г. > р 
` P "2 ба Оо = mH WH OZ? (3E 1 > н > су ti t4 «13 H > 
Ë + U ы LES А 04 «а ө а” H н £4 нн f £4 fa tuu t4 
м |2222 OQ FEAT мч B) Ett» Bap) E: o 1 Е Gz + O 
x «ою ра д, ДО) 0702 ың UH м ө) m н ж се: >; e n kA Z 
u Ë Ë ч О уг ончо зс moH HW ос: м н its г) tr} m ғ“ H 
Ел (ыы ОО "папам Au VHH 2 [o £4 m. ш о F4 £4 
"Ax U ао, е”, Met (прац OEA 30 Су. іл F4 о с D) or 5 N ин 
ае Net wb ыш ин mH AN H H эн E т бы H H 
Hite ош батасы) Осы H >J tn п, B Ë) £4 gi fu 
X ow [x D Dm HMUDDAmOAN DA E4 س‎ н он m uni E! 
v) > ео Ыы ов Д І > < kE > 2) «с ей > P. ul FAT = 
KIRI UZ mn emaudaabnHomomudHze O mm на m m З uj > m о 
о 4 Una et A UT ZU emn к о ta ex о «4 е» га 
9 > t5 und EA п; Сс чо ыт © HH фо OH o4 t >, «О £4 о 
чех = н “СС BAUGQOYU som ` NH WZ Въ <f > п T+ > > 
= на че тыш: гада # +4 piz DH zig Зас Е+ 01 C [x f] H 
22 E44 U HIG OEE! SOMO KS HO а о 400 uU z Г tá 
3 WOR => OOm'doourrmtu жа а; UH Om HO 44 оно н> H wr о ш 
B HANA OF HHOOMEZMOO s # С) fae HO ыға ° Баға r > + Ee O 
$ 2 + а ED Усон панч О сас (9 8 (Е xU) ә о Е он Ез n O tw 
X & NUH аза) ц" Bl. 131 £4 214 гч наб? г Fit ч 
£ < Supr ac Ж (кожа и uU BH DpZ p NDL шы О ыса р 
ш Wh ае чч < ломом сж e Olt n сц каг он “ш 
Фо ыы жи Ех # “Ж ш с са 1 ш о рч fu) = < 
< - iet OOOOQOcQOnuoouoOootupme >t (1 > Ез Ut: => 22 
HORM (nnmrmagnmuánnm mm PHD есу) эц до е «uU m но «ў OH SS 3 
KS Z tí дц. ц М. 3 мы HO "ini ош OO O іра b) ©су] шо oO 
IZH о инн: Шә] бї} су > н Pita жі > H> но N tata r> m 
м Q Юр DLD NTO o ал.» 
зе сэ самого оо «С M А Or O> ОО A ож О eu зо ао со A 
оз CO ОООООСО от ся EAA мэт AH 2A A weed име) oT MID HDA A 
с, сс? >з а сес сес toO се Cc 


ACE AREA OR OLE 
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ER THIS FITTING (Y/N)? 


Т 
>> YOU ARE WORKING CN FITTING UNSEEN о 


RANT TO USE THE DD963 ТҮРЕ 
(Y/N) ? 


Y CONDITION 


S74 5 
DO YOU 
IHE DR 


іт} 


DD STEAL TO DUCTS TI ЧАТ ВЕ ROUND 


ON, IS TNE DUCT CIRCULARS OF RECTANGULAR 


ту 


eee Е ЕВЕ С5055-5ЕСТІОНАУ, 


Ow (а 


ND LIMENSION (FEET) 


ER TEE: LENGTH OF THIS DUCSECILON. (FEET) 


-4 tu: (Оба ше “Уз 
Un 


SG {OU WANT IQ SN MER IHIS FITIING (Y/X)?*? 
»» YCU ARE WORKING CN FITTING NUMBER »» 312204 


00 
Eo? YOU ARE WORKING CN FITTING NUMBER >> 323101 


14 

YOU HAVE SELECTED THE MAIN SECTION OF A DIVERGING WY2. 

THE AIR TO THE ENGINE SHOULD BE FLOWING THROJGd,THIS SECTION. 
JUST ONE CUESTION WHAT IS THE CROSS-SECTIONAL AREA OF THE 
WAIN SECTION? THIS SHOULD BE THE ARZA JUST DOWNSTREAM OF THE 
MUNCTICN AKD DIRECTS FLOY TO THE ЗМСІЧЕ. IT ALSO SHOULD ВЕ 
ОТНЕ FIRST FITTING OF THE BRANCH. 

81.375 

„>> ҮС ARE WORKING CN FITTING NUMBER >> 323102 

26 

YOU HAYE SELECTED A MULTI-BAFFLE TYPZ SILENCEE. 

PACH BAFFLE HAS A STREAMLINED SHAPE. IT IS THE TYPE 

USED IN THE INLETS OF THE DD963. 

s**FIRST QUESTION, @HAT IS THE GAP SETAEEN THE BAFFLES? 
0.333 

jHHAT IS THE THICKNESS OF THE BAFFLES? 

0.666 

„ПНАТ 15 ТНЕ ЦЕНСТН СЕ THE BAFFLES? 

9.33 

„ПНАТ IS THE DIMENSICN OF THE BAFFLES PARALLEL IO THE GAP? 
p 75 

,FHÀT IS THE DIMENSION Or THE MAIN DUCT ACROSS THE GA2S? 
10.5 

,LAST QUESTION, HOW MANY GAPS ARE THERE? 

11 


ES TOU RANI TO ENTER THIS FITIING (Y/N)? 
pee tCU ARE WORKING CN FITTING NUMBER >> 323103 
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5 
АТ IS THE GREATER UPSTREAM CROSS-SECTION DD Pie Sone 
HAT IS THE LEAST DCWHNSTREAM CROSS-SECTION DIMENSION? 


S.ION, WHAT IS IHE GREATEHEDONNSTHEA. 
CTION DICE NSION: 


nie 


DO YOU WENT ТО ENTER THES РЕ (И 
">> YOU ARE RORKING CN FITTING NOSSER 22 И 


6 

YOU HAVE SELECTED A MITERED, RECTANGULAR CROSS-SECTION, ELE OME 
**PTEST7 OUZSTION, WHAT IS TH2 HEISd OC TBE siapa 
_ (THE DIMENSION PARALLEL TO THE TUBN AXIS) 
6.67 

WHAT IS THE WIDTH OF THE ELBOW CROSS-SECTION? 
, (THE CIMENSICN IN THE PLANE OF THE ZUBN) 
7.75 

LAST QUESTION, WHAT IS THE ANGLE OF THE ELBOW TURN 

(0 - $0 DEGRESS)? .. 


90 
DO YOU WANI TO ENTER THIS FITTING (EVR) 


SS YCU ARE WORKING CN FITTING NUMBER >> 323105 

23 

YOU BAVE SELZCTED AVSCAEEN GSS ITRUGCT EON 2 THERDU GT 
„**Ё1852 QUESTION, WHAT IS THE DUCT C20SS-SECTIONAL AREA? 
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,LAST QUESTION, WHAT IS THE FREE FLOW AREA OF THE SCREEN? 
20 015 

DO YCU WANT TO ENTER THIS FITTING (Y/N)? 

‘>> YOU ARE WORKING ON FIZTING NOMBER >> 323106 
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YOU HAVE SELSCIZD IHE BRANCH SECTION ОЕ 1 DIVERGENT aa 
THe BODULE COOLING ALR SHOULL BE 3BANCHINGI IO EE PPE P 5 
INLET AND zsLOWING THROUGH THIS SECTION... ІНІЗ ЗЕ ЕТЕ 
las] т и. ОЕ Fairs дама на 
**"DLBSI QUESTION, “HAT IS THE ANGLE SEC VEE INS и 
AXIS ANT TEE BRANCH FLOW AXIS (DEGREE cs) 
n 
90 
WHAT IS THE CROSS-SECTIONAL ARES OF THE СОДЕ ЕСЕГЕ: 
SECTICN?  IdIS iS 4HERE 3OTH ENGINE AIR AND COOLING Ев 
„JUST UPSTREAM OF TAXE DRANC HG 
= 
„LAST CUESTION, WHAT IS THE CROSS-SECTIONAL ARES OF THE ЕЕ 
5.761 
DO YOU WANI TO ENTER THIS FITTING (YZ 1) 
Y 





бо? YCU ARE NORKING CN FIITING NUMBER >> 328002 
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OM HAVE SecusClED SIRALIGHIT DUC... IT 4! Br ROUND 
DP RECTANG JLAR. 

kee FIRST QUESTION, IS THE DUCT CIRCULAR OR RECTANGULAR  (C/R) ? 
с 

.THz DUCI IS CIRCULAR, ENTER THE DIAMETER (FZET) 
2.708 

SENTER THE LENGTH OF THIS DUCT SECTION. (FEET) 

2.5 

90 YOU WANT TO ENTER THIS FITTING (Y/N)? 
">> ҮС АПЕК WORKING CK FLITING NUNS ERY D> >. 324003 
00 


E YOU ARE WORKING CN FITTING NOMSER >> 335101 
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ШЯШИНАМЕ 521 ЕСТЕП ЗТВАТОНТ DOCS. Ii ArT BS ROUND CR RECTANGULAR: 
ESSLPISSTOQUESIION, ІЗ ІНг DUCT CIRCULAR OR RECTANGULAR  (C/R) ? 
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>> YOU ARE WORKING CN FITTING NUMBER >> 335102 


DO YOU WANT TO ENTER THIS FITTING (Y/N)? 
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Bo YCU ARE WORKING CN FITTING NUMBER >> 3585001 
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CONA E SELECTED TEE GAS TURBINE NODULE AS А 2BARI OF 

EH: COOLING FLOW PASSAGE. NOMCUESTIOIS, JUST NEEDED 

во НО» ЯН YOU WANTED IBE SODULE. 

ЕП ТОУ ТАЯП TO ENTER THIS ЕІТТІІЗС (1/4)? 
>> YOU ARE WORKING CN FITTING NUMSER >> 345002 
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I RAVE SELECTED THE BRANCH ЭЕСТТОМ OF A CONVERGENT 

ШІЕ- TEE НОТ NODULE COOLING AIR SHOULO BE OIALNG THE MAIN 
TENE ЕХНАЦӘТ IN THIS AYE. THIS PFIZIZING SEOULD BE THE LAST 
SITING IN THE ERANCH. 

ЕВЕ 9М557208, ШНАТ IS TPE ANGLE BZIAEEN ІНЕ MAIN FLOW 

AXIS AND THE BRANCH AXIS (DEGREES) ? 


..) 


(e 





АНАТ IS THE CBOSS-SICTICHIAL ABRFA OF T4Z COMBSINES FLC4 

S2CTION? “HIS IS AHERE ZINGINE EXAAUSI AND ЧОРУ ТЕ ССО ОГ 

“ALO JUST DOWNSTREAM OF THEERANCH. 

30.46 | 
LAST QUESTION, WHAT IS THE CROSS-SECTIONAL AREA OF THE 

{ERA NCA? 
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DO YOU WANT ТО ЕЧТЕР ТНІ5 ЕРІТІЛС (ҮІ, | 
E YCU ARE WORKING CN FITTING NUMBER >> 345003 
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„>> YOU ARE WORKING CN FITTING NUMBER >> 356221 

. 1 а 
YOU BAVE SZLECTED À CIRCULAE CONTRACTION. 

,**PIRST QUESTION, WHAT IS THE LENGTH OP THE CONTRACTION? 
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.SHiT IS THRZ UPSTREAN DIAXETER? 

6.237185 

THAT IS THE DOWNSTREAM DIA METER? 

5.4667 

20 YCU WANT IO ENTER THIS FITTING (Y/N)? 
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YOU HAVE SELECTED STRAIGHT DUCT. ІТ ЧАЎ ВЕ ROUND Ca RECTANGULAR. 
+x FIRST QUESTION, IS THE DUCI CIRCULAR Ok RECTANGONTAD SEE 
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БЫ, ЕГО т TRAE СЕ. 5-СТОНАЬ ОТА МЕТЕН ? 
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U SUES ION, АВЕ ӨРТІЛІН NUMBER OF CONCENTRIC VANES 
TALLED TO REDUCE RESISTANCE AND TURBULANCE Ч) > 


DOR OJ PANIIT O ENTER fais FETIING (Y/N) 2 
>> YOU ARE WORKING ON FITTING NUMBER >> 356206 
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fou HAVE SELSCTED STEATGHI DUCT. IT SAN SE ROUND OR RECTANGULAR. 
SORSI ОПЕЗІТОН, ІЗ ІНЕ DUCT CIRCULAR OR RECTANGULAR (C/R) ? 
ВРЕТ 1S с ЕСИ АН ENTER THE DIAMETER (FEET) j 
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OU RANT TO ENTER THIS FITTING (Y/N)? 
>> YOU ARE WORKING CN FITTING NUMBER »» 356207 
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YOO HAVE SELECTED STRAIGHT DUCT. IT MAY BE ROUND OR RECTANGOLAR. 
See ТВТ QUESTION, I5 THES DUCT CIECULAR Ox RECTANGULAR  (C/R) ? 
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DO YCO WANT TO ENTER THiS FITTING (Y/N)? 


Е 
>> YOU ARE WORKING CN РЕШЕНИЕ osu 


THE CONTRACTION? 
I IS TRE UPSTREAM DIAMETER? 
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T IS THE DOWNSTREAM DIAMETER? 


YOU WANT TO ENTER THIS FITTING (Z) 
TCU ARE WƏRX ING CN FIDIINGOUUMHER 25293550200 


U AANI DO ENTEE THAIS FEITTING OLN? 
U ARE WORKING CNX FITTING NUIDER >> 35521] 


O Om UT 


“НАТ SE 
TOU MAY 
S1090? 
DO YCU WANT TO COMPUTE WITH TBE FILE OR QUIT (CZ)? 
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Ee EDITING THE DUCT DATA FILE 


This section demonstrates the  editinj capability of the 


program. Ihe  editcr will Бе demonstrated bby changing a 
ЕЕВС. The fitting chosen 15 ап elbow in the exhaust 
duct. It has cascaded turning vanes installed. By using 


the editor the turning vanes will be removed and an ordinary 
nitered round elbow will be substituted. Any fitting that 


also serves the purpose could be substituted as well. 


The program can also add or delete a fitting. It is 
somewhat limited in the addition  arility. The program can 
not add a fitting to the first of a branch in one step. То 


add a fitting to the duct data file select the index of the 
fitting in the file that the fitting is to be placed after. 
The program will ask what fitting is to be added and then 
ШЕ ч<ес сад enter the fitting directly or from the menu. 
То ааа а fitting at the first of a branch, first add the 
same first fitting presently in the branch after itself, 
then chanye the same index fitting as the first step to the 
desired new first fitting. 

It should be emphasized that the editor does not change 
a system class. If the user wants a different duct arrange- 


ment a new file will have to be entered. 
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DO YCU 


LIV PORC GRISE ANOTHER FITTING (Y/N) ? 
n 
RIT TO NACE АН ОВО ТВЕЕ СНА ЧСЕЗ (Y/N)? 


n 
БИН ЗЕЕ НЫЕ ЖООГО TOU СТАЕ TO GIVE THIS DUCT DATA TILE? 
,YOU MAY USE UP TO A SIX DIGIT INTEGER NUMBER. 
510902 

DO YOU WANT TO COMPUTE WITH THE FILE OR QUIT (C/Q)? 
q 


F. COMPUTING SYSTEM PERFORMANCE 


This section also contains a recorded terminal sess sS 


The ccmputing secticn of the program was exercised here. 
The session has been annotated to point out program 


features. 
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GLOBAL TXTLIB CMSLIB FORTMOD2 MOD2EZH IMSLSP WONIMSL 
LOAD THESIS ( START 
EmSGUTION SEGINS... | 
A ONE-DIMENSIONAL MCDEL FOR THE SYSTEM PERFORMANCE 
OF A MARINE GAS TURBINE INSTALLATION 
BY LCDR. STEPHEN 4. ZZZELL 
VERSION 1.C (“ARCH 30, 1984 
OPTICNS: EUILD À DATA FILE aEPRZSENTING THE DUCT SYSTEM 
EDIT OR CHANGE THE DUCT DATA FILE 
COMPUTE SYSTEM PERFORMANCE 
METHOD: INTERACTIVE INPUT OF DATA, BRANCHING TO DESIRED 
CPTION BY ANSWERING QUzSTfONS 
*#** WARNING, THO NOLL ENTRIES ON NUMERICAL INPOT WILL *** 
*** KILL TdE PROGRAM. ** 
PIRST QUESTION: 
DO YCO HAVE A DATA FILE OF DUCT FITTINGS (Y/N)? 


SO TOUTE TO EDIT THE FICE OR USE IT POR COMPUTATION (E/C)? 
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PORTICN OF IHE EROGRAB INPUTS IHE ENVIROMENTAL CONDITIONS. 
TS aoe IEJTPSRAXWEURE (DEGREES Р)? 
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HESDEPFAWUED SPECEICATIONS. (Y/N)? 
NG YOU DESIRE. 
| ER? 
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HO FU 


ВЕЕТ РОЯ ЕЯ ТОВВЕНЕ SPEED (RPM)? 
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TM ENTERED 
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DE 
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ж с? 


AE TS POVER IJURBINE SPEED (RPM)? 


Кого 
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RESULTS OF TRIS RUN HAVE SEEN ENTERED 
OS FUELLE CALLED "CUIPUI БАТА". 
ВОО АНТ TO COMPUTE WITH DIFFERENT OPERATING CONDITIONS (Y/N)? 
т 


n 
uur TOU RANT TO EDIT TITHE DUCT DATA FILS OR QUIT (Е/0)? 
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с. EXAAINING THE OUTPOT 


Included in tnis section are copies of two files. 
first isa copy of the file the author |) bungee 
Arleigh Burke class example. Тһе other one is 4a copy of 
results from the runs made in the compute section using 


sample rile at two operating points. 
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